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1 Introduction 

Most stars form in dense star clusters deeply embedded in residual gas. The 
populations of these objects range from small groups of stars with TV « dozen 
binaries within a volume with a typical radius of r w 0.3 pc through to 
objects formed in extreme star bursts containing TV « 10^ stars within 
r « a few dozen pc. Star clusters, or more generally dense stellar systems, must 
therefore be seen as the fundamental building blocks of galaxies, whereby a 
differentiation of the term star cluster from a spheroidal dwarf galaxy becomes 
blurred near N ^ 10^ Mq. Both are mostly pressure-supported, that is, ran- 
dom stellar motions dominate any bulk streaming motions such as rotation. 
The physical processes that drive the formation, evolution and dissolution 
of star clusters have a deep impact on the appearance of galaxies. This im- 
pact has many shades, ranging from the properties of stellar populations such 
as the binary fraction and the number of type la and type II supernovae, 
through the velocity structure in galactic disks such as the age-velocity dis- 
persion relation, through to the existence of stellar halos about galaxies, tidal 
streams and the survival and properties of tidal dwarf galaxies, the existence 
of which challenge current cosmological perspectives. Apart from this cosmo- 
logical relevance, dense stellar systems provide unique laboratories in which 
to test stellar evolution theory, gravitational dynamics, the interplay between 
stellar evolution and dynamical processes, and the physics of stellar birth and 
stellar feedback processes during formation. 

Star clusters and other pressure-supported stellar systems on the sky 
merely offer snap shots from which we can glean incomplete information. 
Because there is no analytical solution to the equations of motion for more 
than two stars, these differential equations need to be integrated numerically. 
Thus, in order to proceed to an understanding of these objects in terms of 
the above issues, the researcher needs to resort to numerical experiments in 
order to test various hypotheses as to the possible physical initial conditions 
(to test star-formation theory) or the outcome (to quantify stellar populations 
in galaxies, for example). Initialising a pressure-supported stellar system such 
that the initial object is relevant for the real physical universe is therefore a 
problem of some fundamental importance. 

In the following, empirical constrains on the initial conditions of star clus- 
ters are discussed, and some problems star clusters are relevant for are raised. 
Section [2] contains information on setting-up a realistic computer model of a 
star cluster, including models of embedded clusters. The initial mass distri- 
bution of stars is discussed in section [3l and section |4] delves on the initial 
distribution functions of multiple stars. A brief summary is provided in sec- 
tion [S] 
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1.1 Embedded clusters 

In this section an outline of some astrophysical aspects of dense stellar systems 
is given, in order to help differentiate probable evolutionary effects from initial 
conditions. A simple example clarifies the meaning of this: An observer may 
see two young populations with comparable ages (to within one Myr, say). 
They have similar observed masses but different sizes, a somewhat different 
stellar content and different binary fractions. Do they signify two different 
initial conditions derived from star formation, or can both be traced back to 
a t = configuration which is the same? 



Preliminaries 

Assume we observe a very young population of N stars with an age Tago, and 
that we have a rough estimate of its half-mass radius, ro.5, and embedded 
stellar mass, Meci0 The average mass is 

Also assume we can estimate the star-formation efficiency (SFE), e, within a 
few To 5 for this object, 

e= — , (2) 

Afccl + Mgas 

where Mgas is the gas left over from the star-formation process. The tidal 
radi us of the embedded cl uster can be estimated from the Jacobi limit (eq. 7- 
84 in Binnev &: Tremainel[l987i) as determined by the host galaxy and ignoring 



any contributions by surrounding molecular clouds, 

V 3Mgai y ' ^ ' 

where Mgai is the mass of the Galaxy within the galacto-centric distance D of 
the cluster. This radius is a rough estimate of that distance from the cluster at 
which stellar motions begin to be significantly influenced by the host galaxy. 

The following quantities that allow us to judge the formal dynamical state 
of the system: the formal crossing time of the stars through the object can be 
defined as 

tc.-^, (4) 
a 

wher^ 



^ Throughout this text all masses, m, M, etc. are in units of Mq, unless noted 
otherwise. "Embedded stellar mass" refers to the mass in stars at the theoretical 
time before residual gas expulsion but when star formation has ceased. 

^ As an aside, note that G = 0.0045 pc^/Mq Myr^ and that Ikm/s = 1.02pc/Myr. 
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(5) 

is, up to a factor of order unity, the three-dimensional velocity dispersion of 
the stars in the embedded cluster. Note that these equations serve to estimate 
the possible amount of mixing of the population: If Tago < tcr then the object 
cannot be mixed and we are seeing it close to its initial morphology. It takes a 
few for a dynamical system out of dynamical equilibrium to virialise back 
to dynamical equilibrium. This is not to be mistaken for a relaxation process. 

Once the stars orbit within the object, they exchange orbital energy 
through weak gravitational encounters and rare strong encounters, and the 
system evolves towards a state of energy equipartition. T he energy equ iparti- 



tion time-scale, tms, between massive and average stars (jSpitzen 119871 p. 74), 
which is an estimate of the time massive stars need to sink to the centre of 
the system through dynamical friction on the lighter stars, is 

TO 

^ms — ^rclax- (6) 

^max 

Here, Wmax is the massiv e-star mass, and the chara cteristic two-body relax- 



ation time (e.g. eq. 4-9 in lBinnev fc Tremaindll987l ) is 



N 

ticUx = 0.1 - — — ■ tcr. (7) 

InA' 

This formula refers to a pure TV— body system without embedding gas. A 
ro ugh estimate of trsiax Brn h for an embedded cluster can be found in eq. 8 
of Adams fc Mverd ( 2001 ). The above eq. [7] is a measure for the time a star 
needs to change its orbit significantly away from the initial trajectory and is 
often estimated by calculating the amount of time that is required in order to 
change the velocity of a star, v, by an amount Av « v. 

Thus, if for example, Tago > icr and Tago < irciax, then the system is 
probably mixed and close to dynamical equilibrium, but it is not relaxed yet, 
i.e., did not have sufficient time for the stars to exchange significant amount 
of orbital energy. Such a cluster may have erased its sub-structures. 

Fragmentation and size 

The very early stages of cluster evolution on a scale of a few pc are dominated 
by gravitational fragmentation of a turbulent m agnetize d contract ing molec- 
ular cloud core ( Clark e, Bon nell fc Hillenbrand 2000; Mac Low fc K lcssenl 



120041: iTillev fc Pudritz Il2007i) . The existing gas-dynamical simulations show 



the formation of contracting filaments which fragment into denser cloud cores 
that form sub-clusters of accreting proto-stars. As soon as the proto-stars emit 
radiation and outflows of sufficient energy and momentum to affect the cloud 
core these computations become expensive as radiative transport and deposi- 
tion of momentum and mechanical energy by non-isotrop ic outflows are dif- 



ficult to handle with given present co mputational means (jStamatellos et al. 



20071 : iDale. Ercolano fc Clarke 1120071 ). 
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Observations of the very early stages at times ^ few 10^ yr suggest proto- 
clusters to have a hierarchical proto-stellar distribution: a number of sub- 
clusters with radii ^0.2pc and separated in velo city space are oft en seen 



embedded within a region less than a pc across (|Testi et al. II2OOOI ). Most 



of these sub-clusters may merge to form a more massive embedded cluster 
(IScallv fc Clarke II2002I: [Pellhauer fc Kroupa II2005I) . 

It is unclear though if sub-clumps typically merge before residual gas blow- 
out or if the residual gas is removed before the sub-clumps can interact sig- 
nificantly, nor is it clear if there is a systematic mass dependence of any such 
possible behaviour. 



Mass segregation 

Whether or not star clusters or sub-clusters form mass-segregated remains an 
open issue. Mass segregation by birth is a natural expectation because proto- 
stars near the density maximum of the cluster have more material to accrete. 
For these, the ambient gas is at a higher pressure allowing proto-stars to ac- 
crete longer before feedback termination stops further su bstantial gas inflow 



and t he coagulation of pr otostars is more likely there (jZinnecker fc Yorke 



2007t iBonnell et al.l 120071 ). Ini tially mass-segregated sub-clusters preserve 



mass segregation upon merging ( McMillan. Versperini fc Portegies Zwart I2OO7I 1 . 
However, for m/mmax = 0.5/100 and N^5 x 10'^ stars it follows from eq. [H] 
that 

i.e., a 100 Mq star sinks to the cluster centre within roughly a crossing time 
(see table [T] below for typical values of tcr)- 

Currently we cannot say conc lusively if mass segregation is a birth phe- 



nomenon (e.g. ICouliermis et al.~ 2004), or whether the more massive stars 



form anywhere throughout the proto-cluster volume. Star clusters that have 
already blown out their gas at ages of one to a few Myr are typically mass- 
segregated (e.g. R136, Orion Nebula Cluster). 

Affirming natal mass segregation would impact positively on the notion 
that massive stars {^10 Mq) only form in rich clusters, and negatively on 
the sugg estion that they can also form in isolation (for recent w ork on this 
topic see lLi. Klessen fc Mac Low II20Q3I: IParker fc Goodwin Il2007^ . 



Feedback termination 



The observationally estimated SFE (eq. [2]) is (jLada fc Lada 1 120031 ) 

0.2<e<0.4 (9) 

implying that the physics dominating the star-formation process on scales less 
than a few pc is stellar feedback. Within this volume, the pre-cluster cloud 
core contracts under self gravity thereby forming stars ever more vigorously, 
until feedback energy suffices to halt the process {feedback-termination) . 
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Dynamical state at feedback termination 

Each prot o-star needs about t^s ^ 10^ yr to accumulate about 95 per cent of 
its mass ( Wuchterl fc Tscharnuter 20031 ). The proto-stars form throughout 



the pre-chister volume as the proto-cluster cloud core contracts. The overall 
pre-cluster cloud-core contraction until feedback-termination takes (eqs [H [5]) 

2 /MeclV^ I 
ici.form ~ few X -j= [ rjg, (10) 



G 

(a few times the crossing time), which is about the time over which the cluster 
forms. Once a proto-star condenses out of the hydro-dynamical flow it becomes 
a ballistic particle moving in the time-evolving cluster potential. Because many 
generations of proto-stars can form over the cluster-formation time-scale and 
if the crossing time through the cluster is a few times shorter than ici,formi 
then the very young cluster is mostly virialised when star formation stops and 
at the removal of the residual gasH It is noteworthy that for ro.s = 1 pc 

ips^ici.form for ^ > lO^'^ Mq (11) 

(the proto-star formation time formally surpasses the cluster formation time) 
which is near the turnover mass in the old-star-cluster mass function (e.g. 
lBaumgardt1ll998f ). 



A critical parameter is thus the ratio 

^cl.form /-, (..x 

r = — . (12) 

If it is less than unity then proto-stars condense from the gas and cannot 
virialise in the potential before the residual gas is removed. Such embedded 
clusters may be kinematically cold if the pre-cluster cloud core was contract- 
ing, or hot if the pre-cluster cloud core was pressure confined, because the 
young stars do not feel the gas pressure. 

In those cases where r > 1 the embedded cluster is approximately in virial 
equilibrium because generations of proto-stars that drop out of the hydro- 
dynamic flow have time to orbit the potential. The pre-gas-expulsion stellar 
velocity dispersion in the embedded cluster (eq. [5]) may reach a — 40pc/Myr 
if A/eci — 10^'^ M0 which is the case for ero.5 < Ipc. This is easily achieved 
since the radius of one-Myr old clusters is ro.5 ~ 0.8 pc with no dependence o n 
mass. Some observationally explored ca ses are discussed by Krouoa I ( 20051 ). 



Notably, using K-band number counts, Gutermuth et al. (|2005l ) appear to 



find evidence for expansion after gas removal. 



^ A brief transition time ttr <S tci.form exists during which the star-formation rate 
decreases in the cluster while the gas is being blown out, but for the purpose of 
the present discussion this time may be neglected. 
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Interestingly, recent Spitzer results suggest a scaling of the characteristic 
projected radius R with massQ 

Mecl cc (13) 



(| Allen et al. Il2007f ). so the question how compact embedded clusters form 
and whether there is a mass-radius relation needs further clarification. Note 
though that such a scaling is obtained for a stellar population that expands 
freely with a velocity given by the velocity dispersion in the embedded cluster 
(eq.©, 

r{t) f^r^ + at A^.i = 1 (^^^^^) , (14) 

where Tq ^ 1 pc is the birth radius of the cluster. Is the observed scaling then a 
result of expansion from a compact birth configuration after gas expulsion? If 
so, then it would require more-massive systems to be dynamically older, which 
is at least qualitatively in-line with the dynamical time-scales decreasing with 
mass. Note also that the observed scaling (eq. [13]) cannot carry through to 
Meci ^ 10'' Mq because the resulting objects would not resemble clusters. 

There are two broad camps suggesting on the one han d side that molec- 
ular clouds and star clusters form on a free-fall time-scale ( Elmegreen 2000l : 



Hartmann 1 12003' 'Elmegreen"2007^ and on the other that many free-fall times 



are needed (jKrumholz fc Tan , ,20071) . The former implies r « 1 while the 



latter implies r > 1. 

Thus, currently unclear issues concerning the initialisation of A^-body mod- 
els of embedded clusters is the ratio r, and whether a mass-radius relation 
exists for embedded clusters before the development of HIT regions. To make 
progress I assume for now that the embedded clusters are in virial equilibrium 
at feedback termination (r > 1) and that they form highly concentrated with 
r ^ 1 pc independently of mass. 

The mass of the most massive star 

Young clusters show a well-defined correlation between the mass of the most 
massive star, mmax, in dependence of the stellar mass of t he embedded cluster 



Meci, whic h appears to saturate at mmax* ~ 150 (jWeidner fc Kroupa 



2004Ll2006h . This IS visualised in fig.[TJ This correlation may indicate feedback 



termination of star formation within the proto-cluster volume coupled to the 
most mass ive stars forming latest, or turning-on at the final stage of cluster 

formation ([Elmegreen I [nil). 



The evidence for a universal upper mass cutoff near 

TOmax*~150MQ (15) 



* Throughout this text, projected radii are denoted by R, while the 3D radius is r. 
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3 4 
log,o(Meci/Me) 



Fig. 1. The maximum stellar mass, 



as a function of the stellar mass of the 



embedded cluster, Meci (Weidner, private com munication: an updated version of the 
data presented in IWeidner fc Kroupa" 2006h. T he solid triangle is an SPH model 
of star-cluster formation bv lBonnell et aL ( 20031 ). while the solid curve stems from 
stating that there is exactly one most massive star in the cluster, 1 = dm 

with the condition Meci = J^^'^^ m ^{m) dm, where is the stellar IMF. The 

solution can only b e obtained numerically but an easy to use wel l-fitting function 
has been derived bv lPflamm-Altenburg. Weidner fc Kroupal (|2007l ). 



(IWeidner k Kroupa I l200i: iFieer 1 120051: lOev fc Clarkel 120051: iKoen I l2006t 



Mai'z Apellaniz et al.l 2007 : Zinnecker fc Yorke 20071 ) seems to be rather 



well established in populations with metallicities ranging from the LMC 
{Z w 0.008) to the super-solar Galactic centre (ZjJO.02) such that the 
stellar mass function (MF) simply stops at t hat mass. This mass need s to 
be understood theoretica lly (see discussion in iKroupa fc Weidner llioosi and 
Zinnecker fc Yorke1l2007t) . It must be a result of stellar structure stability, but 



may be near 80 M© as predicted by theory if the most massive stars re side 
in near-equal component-mass binary systems ( Kroupa fc Weidnei~ll2005fl . It 
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may also be that the calculated ste llar masses are significantly overestimated 



(jMartins. Schaerer fc Hillier Il2005f) 



The cluster core of massive stars 

Irrespectively of whether the massive stars (JJlOAf©) form at the cluster 
centre or whether they segregate there due to energy equipartition (eq. [6]), 
they ultimately form a compact sub-population that is dynamically highly 
unstable. Massive stars are ejected from such cores very efficiently on a core- 
crossing time-scale, and for example the well-studied Orion Nebula cluster 
(ONC) has probably already shot out 70 p er cent of its stars more massive 
than 5Mq (jPflamm-Altenburg fc Kroupa i2006). The prope rties of O and 
B runaway stars have been used by Clarke fc Pringle ( 19921 ) to deduce the 
typical birth configuration of massive stars, finding them to form in binaries 
witli similar- mass components in compact small- groups devoid of low-mass 
stars. Among others, the core of the Orion Nebula Cluster (ONC) is just such 
a system. 



The star-formation history in a cluster 

The detailed star-formation history in a cluster contains information about 
the events that build-up the cluster. Intriguing is the recent evidence for some 
clusters that while the bulk of the stars have ages differe nt by less than a 



few 1 0^ yr, a small f raction of older sta rs are often harboured (jPalla fc Stabler 



2000i for the ONC. ISacco et al.ll2007l for the cr Orionis cluster). This may be 
interpreted to mean that clusters form over about 10 Myr with a final highly 
accelerated phase, in support of the notion that t urbulence of a magnetiz ed 
gas determines the early cloud-contraction phase (jKrumholz fc Tan 1120071 ). 

A different interpretation would be that as a pre-cluster cloud core con- 
tracts on a free-fall time-scale, it traps surrounding field stars which thereby 
become formal cluster members: Most clusters form in regions of a galaxy 
that has seen previous star formation. The velocity dispersion of the previ- 
ous stellar generation, such as an expanding OB association, is usually rather 
low, around a few km/s to 10 km/s. The deepening potential of a newly- 
contracting pre-cluster cloud core will be able to capture some of the preceding 
generation of stars such that these older star s become formal cluster mem- 



bers a lthough they did not form in this cluster. iPfiamm-Altenburg fc Kroupa 



( 20071) study th i s prob lem for the ONC showing that the reported age spread 
bv IPalla et aTl (|2007l ) can be accounted for in this way. This suggests that 
the star-formation history of the ONC ma y in fact not have started about 
10 Myr ago, supporting the argument by El megreen I ([2000); Elmegreen (20o3) 
and iHartmann I (j2003l ) that clusters form on a timescale comparable to the 
crossing time of the pre-cluster cloud core. Additionally, the sample of clus- 
ter stars may be contaminated by enhanced fore- and back-ground densi- 
ties of field stars due to focussing o f stellar orbits during cluster formation 
(jPfiamm-Altenburg fc Kroupa I l2007l ). 
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Fo r very massive clusters such as uj Cen, iFellhauer. Kroupa fc Evans . 
( 2006^ show that the potential is sufficiently deep such that the pre-cluster 
cloud core may capture the field stars of a previously existing dwarf galaxy. 
Up to 30 per cent or more of the stars in oj Cen may be captured field stars. 
This would be able to explain an age spread of a few Gyr in the cluster, and 
is consistent with the notion that u Cen formed in a dwarf galaxy that was 
captured by the Milky Way. The attractive aspect of this scenario is that 
u! Cen need not have been located at the center of the incoming dwarf galaxy 
as a nucleus, but within its disk, because it opens a larger range of allowed 
orbital parameters for the putative dwarf galaxy moving about the Milky 
Way. The currently preferred scenario in which to Cen was the nucleus of the 
dwarf galaxy implies that the galaxy was completely stripped while falling 
into the Mi lky Way leaving only its nucleus on its current retrograde orbit 
( Zhaoll2004 . The new scenario allows the dwarf galaxy to be absorbed into 
the Bulge of the MW with uj Cen being stripped from it on its way in. 

Another possibility for obtaining an age spread of a few Gyr in a massive 
cluster such as uj Cen is gas accretion from a co-moving inter-stellar medium 



(jPflamm-Altenburg fc Kroupa 1120081 ) . This could only have worked for oj Cen 



before it became unbound from its mother galaxy, though. That is, the cluster 
must have spent about 2-3 Gyr in its mother galaxy before it was captured 
by the Milky Way. 

This demonstrates beautifully how an improved understanding of dynam- 
ical processes on scales of a fewpc impinges on problems related to the for- 
mation of galaxies and cosmology (through the sub-structure problem). 

Finally, the increasingly well-documented evidence for stellar populations 
in massive clusters with different metallicities and ages, and in some cases 
even significant He enrichment, may also suggest secondary star formation 
occurring from material that has been pre-enriched from a previous generation 
of stars in the cluste r. Different IM Fs need to be invoked for the populations 
of different ages (see Piotto 20081 for a review). 



Expulsion of residual gas: 

When the most massive stars arc O stars they destroy the proto-cluster neb- 
ula and quench further star formation by first ionising most of it (feedback 
termination). The ionised gas, being now at a temperature near 10^ K and in 
serious over-pressure, pushes out and escapes the confines of the cluster vol- 
ume with the sound speed (near lOkm/s) or faster if the winds being blown 
off O stars with velocities of thousands of km/s impart sufficient momentum. 

There are two analytically tractable regimes of behaviour: instantaneous 
gas removal and slow gas expulsion over many crossing times: 

• Instantaneous gas expulsion, Tgas = 0: The binding energy of the object of 
mass M is 

G Af ^ 1 

^ci,bind = + -Ma2 <0. (16) 

r 2 
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Before gas expulsion, M = Minit — -A/gas + Afeci — > (becomes) M, and 

f^init = — ' (17) 

' init 

After instantaneous gas expulsion, Maftor = — * M, but aaftcr — 
Cinit — * o" and we get the new binding energy 

-£'cl, bind, after = + -^after '^fnif (1^) 

nnit 2 

But the cluster virialises into a new equilibrium, such that, by the scalar 
virial theoren|f| 

1 GMafter .„^ 

^cl, bind, after — —7: , 

^ ''after 

and on equating these two expressions for the final energy and using eq. [17] 
it follows that 

^after ^cc\ /ori\ 

~ Mecl - Mga. ■ 

thus, as Mgas — > Mcci, i-c. e — > 0.5 from above, rafter — > 
oo. This means that as the SFE approaches 50 per cent from above, 
the cluster unbinds. But by eg. [H this result would imply either (see 
Kroupa. Aarseth &: Hurlev 200ll and references therein) 



- all clusters with OB stars (and thus Tgas <^ icr) do not survive gas 
expulsion, or 

- the clusters expel their gas slowly, Tgas 3> icr, which may be the case if 
surviving clusters such as the Pleiades or Hyades formed without OB 
stars. 

Slow gas removal, Tgas ^ ^cr, ''gas — > oo: By eg . and assuming that 
an infinitesimal mass of gas is removed instantaneously, 

^init Mnit - <5 A/gas - SMgs.s 

For infinitesimal steps and taking for convenience dM < but dr > 0, 

r-dr M + dM 



M + 2 dM 



(22) 



Re- arranging. 



so that 



^ = ^^fl-2^^...), (23) 
r M \ M ' ' ^ ' 



^ln!:^.ln^, (24) 

r M rinit A/aftor 



The scalar virial theorem: 2K + W ^ Q E ^ K + W ^ (1/2) W, where 

K, W axe the kinetic and potential energy and E is the total energy of the system. 
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upon integration of the differential equation. Thus, 

?^after _ -^^ccl + -^^gas _ 1 ^25) 

and for example for a SFE of 20 per cent, the cluster expands by a factor 
of five, rafter = 5rinit, without dissolving. 

Table [1] gives an overview of the type of behaviour one might expect for 
clusters with increasing number of stars, TV, and stellar mass, Mgd and for 
two characteristic radii of the embedded stellar distribution, ro.5. It can be 
seen that the gas-evacuation time scale becomes longer than the crossing 
time through the cluster for Mod ^ 10^ Mq. Such clusters would thus undergo 
adiabatic expansion as a result of gas blow out. Less-massive clusters are 
more likely to undergo an evolution that is highly dynamic and that can 
be described as an explosion (the cluster pops). For clusters without O and 
massive B stars, nebula disruption probably occurs on the cluster-formation 
time-scale, « 10^ yr, and the evolution is again adiabatic. A simple calculation 



Mecl/M0 


N 


stars? 


tcr/Myr 




tcr/Myr 








(ro.5 = 


0.5 pc 


0.5 pc 


1 pc 


Ipc) 


40 


100 


N 


0.9 




2.6 




100 


250 


Y/N 


0.6 


0.08 


1.6 


0.2 


500 


1250 


Y 


0.3 


0.2 


0.7 


0.1 


10^ 


2.5 X 10^ 


Y 


0.2 


0.25 


0.5 


0.2 


10^ 


2.5 X 10" 


Y 


0.06 


0.8 


0.2 


0.5 


10^ 


2.5 X 10^ 


Y 


0.02 


2.5 


0.05 


2 


10** 


2.5 X 10^ 


Y 


0.006 


8.3 


0.02 


5 



Table 1. Notes: O stars — "Y" if the maximum stellar mass in the cluster surpasses 
8Mq (fig. HI; the average stellar mass is taken to be m = OAMq in all clusters; a 
star-formation efficiency of e = 0.3 is assumed; the crossing time, tcv, is eq. |4l the 
pre-supernova gas evacuation time-scale is Tgas = r/vth, where Vth = lOkm/s is the 
approximate sound velocity of the ionised g = 0.05 Myr for r = 0.5 pc, while 

Tgas = 0.1 Myr for r = Ipc. 



of the amount of energy deposited by an O star into its surr ounding cluster - 
nebula suggests it to be larger than the nebula binding energy ( Kroupa Il2005f ). 
This, however, only gives at best a rough estimate of the rapidity with which 
gas can be expelled; an inhomogeneous distribution of gas leads to the gas 
removal occurring preferably along channels and asymmetrically, s uch that the 



overa ll gas-excavation process is highly non uniform and variable ([Dale et al. 
20051) . 



The reaction of clusters to gas expulsion is best studied numeric ally with 
iV— b ody codes. T he pioneering such experi ments were perfor m ed bylTutukov 
(flQTS) followed bv lLada et air(|l984l ). and (|Goodwinlll997allbl . llQQSh studied 



Dense Stellar Populations; Initial Conditions 



13 



gas expulsion by supernovae from young globular clusters. Fig. [2] shows the 
evolution of an ONC-type initial cluster with a stellar mass Mod ~ 4000 Mq 
and a canonical IMF (eq. I124p and stellar evolution, a 100 per cent initial 
binary population (section [4. 2p in a solar-neighbourhood tidal field, e = 1/3, 
and spherical gas blow-out on a thermal time-scale (vth = 10 km/s). The figure 




-10 12 
logiot [Myr] 

Fig. 2. The evolution of the 5, 10, 20..., 50 per cent Lagrange radii and the core 
radius {Rc = rc, thick lower curve) of the ONC-type cluster discussed in the text. 
The gas-mass is shown as the dashed line: the cluster spends 0.6 Myr in an embedded 
phase before the gas is blown out on a therm al time-scale. The tidal radiu s (eq. |3} 
is shown by the upper thick solid curve. From iKroupa, Aarseth fc Hurlevl (pOOIl . 

demonstrates that the evolution is far more complex than the simple analyt- 
ical estimates above suggest, and in fact a substantial Pleiades-type cluster 
emerges after loosing about 2/3rd of the initial stellar population (see also 
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p.ll6p. Subsequent theoretical work based on an iterative scheme according to 
which the mass of unbound stars at each radius is removed successively shows 
that the survival of a cluster depends not only on e, Tgas/^cr and rtiH, but also 



on th e detailed shape of the stellar distribution function (jBoilv fc Kroupa 
2003f) . For instantaneous gas removal, e « 0.3 is a lower limit for the SFE 



below which clusters cannot survive rapid gas blow-out. This is significantly 
smaller than the critical value of e = 0.5 below which the stellar system be- 
comes formally unbound (eq. [20|) . However, if clusters form as complexes of 
sub-clusters, each of which pop individuall y, then overall cluste r survival is 
enhanced to even smaller values of e ~ 0.2 ( Fellhauer fc Kroupa I [20 05). 

// clusters pop, and which fraction of stars remain in a post-gas expul- 
sion cluster, depends critically on the ratio between the gas-removal time 
scale and the cluster crossing time. This ratio thus mostly defines which clus- 
ters succumb to infant mortality^ and which clusters merely suffer cluster 
infant weight loss. The well-studied observational cases do indicate that the 
removal of most of the residual gas does occur within a cluster-dynamical time, 
Tgas/^cr ^ 1- Examples noted ( Kroupa 20051 ) are the ONC and R136 in the 
LMC both having significant super-virial velocity dispersions. Other exam- 
ples are the Treasure- Chest cluster and the very young star-bursting clusters 
in the massively-interacting Antennae galaxy which appear to have HII re- 
gions expanding at velocities such that the cluster volume may be evacuated 
within a cluster dynamical time. However, improved empirical constraints are 
needed to further develop an understanding of cluster survival. Such observa- 
tions would best be the velocities of stars in very young star clusters, as they 
should sho w a radially expanding stella r population. 

Indeed, iBastian fc Goodwin note that many young clusters have 

a radial-density profile signature expected if they are expanding rapidly, sup- 
porting the notion of fast gas blow out. For example, the 0.5 — 2Myr old 
ONC, which is kn own to be super-virial with a virial mass about twice 
the observed mass ( Hillenbrand fc Hartmannlll998 ). has already expelled its 
resid ual gas and is expanding rapidly t hereby probably having lost its outer 
stars ( Kroupa. Aarseth fc Hurlev 200ll ). The super-virial state of young clus- 
ters makes measurements of their mass-to-light ratia a bad measure of the 
stellar mass within them ( Goodwin fc Bastian 2006[ ). and rapid dynami- 
cal ma ss-segregation likewise makes naive measurements of the M/L ratia 
wrong jBoilv et al.ll2005l: iFleck et al.l l2006f). ICoodwin fc Bastian I (|2006l ) and 
de Griis fc Parmentier ( 2007t ) find the dynamical mass-to-light ratia of young 



clusters to be too large strongly implying they arc in the process of expanding 
aft er gas expulsion . 



Weidner et al. I (|2007t ) attempted to measure infant weight loss by using 



a sample of young but exposed Galactic clusters and applying the maximal- 
star-mass vs cluster mass relation from above to estimate the birth mass of 
these clusters. The uncertainties are large, but the data firmly suggest that 
the typical cluster looses at least about 50 per cent of its stars. 



Dense Stellar Populations: Initial Conditions 



15 



Binary stars 



Most stars form as binaries with, as far as can be stated today, universal orbital 
distribution functions (section |4]). Once a binary system is born in a dense 
environment, it is perturbed which changes its eccentricity and semi-major 
axis, or it undergoes a relatively strong encounter which disrupts the binary 
or hardens it perhaps with exchanged companions. The initial binary popula- 
tion therefore evolves on a cluster crossing time scale, and most soft binaries 
are disrupted. It has been shown that the properties of the Galactic field bi- 
nary population can be explained in terms of the binary properties observed 
for very young populations if t hese go throug h a dense cluster enviroment 
{dynamical population svnthesis \Kvouv8i\l995ai) . A dense cluster environment 
hardens existing bina ries fp|65p which increa ses the SNIa rate in a galaxy with 
many dense clusters ( Shara fc Hurley 2002f) . 

Binaries have been realised to be significant energy sources (see also sec- 
tion 14. ip : a hard binary that interacts via a resonance with a cluster field 
star ejects one star with a velocity Vcj ^ cr. The ejected star either leaves the 
cluster causing cluster expansion such that a drops, or it shares some of its 
kinetic energy with the other cluster field stars through gravitational encoun- 
ters causing cluster e xpansion. Binaries in a cluster core can thu s halt and 
reverse core collapse (jMevlan &: HeggidllQQTt iHeggie fc Hut 1120031 ). 



Mass loss from evolving stars 



An old globular cluster with a turn-off mass near 0.8 Mq will have lost 30 per 
cent of the mass that remained in i t after gas expulsion due to stellar evo- 
lution (|Baumgardt fc Makino]|2003[) . As the mass loss is most rapid during 



the earliest times after re-virialisation after gas expulsion, the cluster expands 
further during this time. This is nicely seen in the Lagrange radii of realistic 
cluster- formation models ( Kroupa. Aarseth fc Hurlev 200lf l. 



1.2 Some implications for the astrophysics of galaxies 



In general, the above have a multitude of implications for galactic and stellar 
astrophysics: 



1. The heaviest-star- 



star-cluster-mass correl ation constrains feedback mod- 
1983f ). It also implies that by 



els of star cluster formation (jElmegreen 



adding up all IMFs in all young clusters in a galaxy, the integrated galax- 
ial initial mass function (IGIMF) is steeper than the invariant stellar 
IMF observed in star clusters with important implications for the m ass- 
metallicity relation of galaxies (jKoeppen. Weidner fc Kroupa 1120071) . Ad- 
ditionally, star-formation rates (SFRs) of dwarf galaxies can be underesti- 
mated by up to three orders of magnitude because Ha-dark star formation 
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becomes possible (|Pflainni-Altenburg. Weidner fc Kroupall2007( ). This in- 
deed constitutes an important example of how sub-pc processes influence 
the physics on cosmological scales. 

The deduction that type II clusters probably pop (p. [TU]) implies that 
young clusters will appear to an observer to be su per-viriah i.e. to have 



a dynamical mass larger than the luminous mass ([Bastian fc Goodwin 



20061 : Ide Grijs fc Parmentier Il2007^ 



3. It also implies that galactic fields can be heated, and may also lead 



to gal actic thick-disks and stellar halos around dwarf galaxies (|Kroupa 



2002bD . 



4. The variation of the gas expulsion time-scale among clusters of differ- 
ent type implies that the star-cluster ma ss function (CMF) is re -shaped 
rapidly, on a time-scale of a few ten Myr (jKroupa fc Boily II2002I ). 

5. Associated with this re-shaping of the CMF is the natural production 
of pop ulation II st ellar halos during cosmologically e arly sta r-formation 



20081 ) 



bursts (jKroupa fc Boily 2002; Parmentier fc Gilmore 2007; B aumgardt. Kroupa fc Parmentier 



The properties of the binary-star population observed in Galactic fields 
are shaped by dynamical encounters in star clusters before the stars leave 
their cluster (section 2]) . 



Points [1H5] are considered in more detail in the rest of section [T] 



Stellar associations, open clusters and moving groups 

As one of the important implications of point [2l a cluster in the age range 
1 — 50 Myr will have an unphysical M/L ratio because it is out of dynamical 



equil i brium rather than having an a bnormal stellar IMF ijBastian fc Goodwin 



equii i prmm ratner tnan navmg an a t 
200d Ide Griis fc Parmentier II2007I ) 



Another implication is that a Pleiades-like open cluster would have been 
born in a very dense ONC-type configuration and that, as it evolves, a 
moving- group- 1 is established during the first few dozen Myr which comprises 
roughly 2 /3rd of the initial stellar population and is expanding outwards with 
a veloc ity dispersion which is a function of the pre-gas-expulsion configu- 
ration teroupa. Aarseth fc Hurley" 2001 ). These computations were in fact 



the first to demonstrate, using high-precision iV-body modelling, that the 
re-distribution of energy within the cluster during the embedded phase and 
during the expansion phase leads to the formation of a substantial remnant 
cluster despite the inclusion of all physical effects that are disadvantageous 
for this to happen (explosive gas expulsion, low SFE e — 0.33, Galactic tidal 
field and mass loss from stellar evolution and an initial binary-star fraction 
of 100 per cent, see fig. [2]). Thus, expanding OB associations may be related 
to star-cluster birth, and many OB associatio ns ought to have remnant star 
clusters as nuclei (see also lClarke et al. 20051 ). 
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As the cluster expands becoming part of an OB association, the radiation 
from its massive stars produce exp anding HII regions that may trigger further 



star formation in the vicinity (e.g. iGouHermis et al. 1120071) 



A moving-group-II establishes later - the classical moving group made-up 
of stars which slowly diffuse/evaporate out of the re-virialised cluster remnant 
with relative kinetic energy close to zero. The velocity dispersion of moving 
group I is thus comparable to the pre-gas-expulsion velocity dispersion of the 
cluster, while moving group II has a velocity dispersion close to zero. 

The velocity dispersion of galactic-field populations and galactic 
thick disks 

Thus, the moving-group-I would be populated by stars that carry the initial 
kinematical state of the birth configuration into the field of a galaxy. Each 
generation of star clusters would, according to this picture, produce overlap- 
ping moving-groups-I (and II), and the overall velocity dispersion of the new 
field population can be estimated by adding in quadrature all expanding pop- 
ulations. This involves an integral over the embedded-cluster mass function, 
^oci(Afcci), which describes the distribution of the stellar mass content of clus- 
ters when they are born. Because the embedded cluster mass function is known 



to be a powe r-law, this integral can be calculated for a first estimate (jKroupa 
2002bl I2OO5I) . The result is that for reasonable upper cluster mass limits in 
the integral, Mcci^lO^M©, the observed age-velocity dispersion relation of 
Galactic field stars can be re-produced. 

This theory can thus explain the much debated energy deficit: namely 
that the observed kinematical heating of field stars with age could not, until 
now, be explained by the diffusion of orbits in the Galactic disk as a resul t 
of scattering on molecular clouds, spiral arms and the bar (jJenkins I ll992h . 



Because the velocity-dispersion for Galactic-field stars increases with stellar 
age, this notion can also be used to map the star-formation history of the 
Milky- Way disk by resorting to the observed correlation between the star- 
formation rate in a galaxy a nd the maximum star-cluster mas s born in the 



population of young clusters (jWeidner. Kroupa fc Larsen 1120041) 



An interesting possibility emerges concerning the origin of thick disks. If 
the star formation rate was sufficiently high about 11 Gyr ago, then star clus- 
ters in the disk with masses up to 10^'^ Mq would have been born. If they 
poppe d a thick disk wi th a velocity dispersion near 40 km/s would result nat- 
urally ( Kroupa 2002bl ). This notion for the o rigin of thick disks appears to be 



qualit atively supported by the observations of lElmegreen. Elmegreen fc Sheets 



( 20041 ) who find galactic disks at a redshift between 0.5 and 2 to show massive 



star-forming clumps. 

Structuring the initial cluster mass function 

Another potentially important implication from this theory of the evolution of 
young clusters is that if the gas-expulsion-time-to-crossing-time ratio and/or 
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the SFE varies with initial (embedded) cluster mass, then an initially feature- 
less power-law mass function of embedded clusters will rapidly evolve to one 
with peaks, dips and turnovers at cluster masses that characterize changes in 
the broad physics i nvolved . 

As an example, [Adami diooO) and lKroupa fc Boilv I (|2002l ) assumed that 



the function 

A/ici = fst A/cci (26) 
exists, where Med is as above. Mid is the classical initial cluster mass and 

fst = /.t(Mcci). (27) 



According to Kroupa fc Boilv I ( 20021 ). the classical initial cluster mass is 



that mass which is inferred by classical A'^-body computations without gas 
expulsion (i.e. in effect assuming e = 1, which is however, unphysical). 
Thus, for example, for the Pleiades, Md « 1000 A/q at the present time 
(age about 100 Myr). A classical initial model would place the initial clus- 
ter mass near Mic\ ~ 1500 Mq by using standard A^-body calculations to 
quantify the secular evaporati on of stars from an initially bound and viri- 



alised classical cluster (Portc gies Zwart et al. 1120011 ). If, however, the SFE was 



33 per cent and the gas-expulsion time-scale was comparable to or shorter 
than the cluster dynamical time, then the Pleiades would have been born 
in a compact configuration resem bling the ONC and with a mass of em- 
bedded stars of Med ~ 4000 Mq (|Kroupa. Aarseth fc Hurlev I l200ll ). Thus, 
/st(4000 Mq) = 0.38 (= 1500/4000). 

By postulating that there exist three basic types of embedded clusters 



S that 
2002h . namely 



Type I: clusters without O stars (Mod ^ 10^-^ M©, e.g. Taurus- Auriga pre- 

main sequence stellar groups, p Oph), 
Type II: clusters with a few O stars (lO^-^ < A/oci/Mq < 10^■^ e.g. the ONC), 
Type III: clusters with many O stars and with a velocity dispersion compara- 
ble to or higher than the sound velocity of ionized gas (Med ~ 10^'^ Mq), 

it can be argued that f^t ~ 0.5 for type I, fst < 0.5 for type II and fst ~ 0.5 
for type III. The reason for the high fst values for types I and III is that 
gas expulsion from these clusters may be longer than the cluster dynamical 
time because there is no sufficient ionizing radiation for type I clusters, or 
the potential well is too deep for the ionized gas to leave (type III clusters). 
The evolution is therefore adiabatic (eq. [25] above). Type II cl usters undergo 



a disruptive evolution and witness a high infant mortality rate (jLada fc Lada I 
[2003), therewith being the pre-cursors of OB associations and Galactic clus- 
ters. This broad categorisation has easy-to-understand implications for the 
star-cluster mass function. 

Under these conditions and an assumed functional form for fst — /st(Afoci), 
the power-law embedded cluster mass function transforms into a cluster 
mass function with a turnover near 10^ A/q and a sharp peak near 10'^ Mq 
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( Kroum k Boilv Il2002f) . This form is strongly remin iscent of the initial glob 



ular cluster mass funct i on wh i ch is inferred by e. g. IVesperini I (jl998l . l200lh : 



Parmentier fc Gilmore ( 2005 ): Baumgardt I ( 1998 ) to be required for a match 
with the evolved cluster mass function that is seen to have a universal turnover 
near 10^ Mq. By the reasoning given above, this "initial" CMF is, however, 
unphysical, being a power-law instead. 

This analytical formulation of the problem has been verified nicely using 
A^-bod y simulations combined with a realist i c trea- tment of residual gas expul- 
sion bv lBaumgardt. Kroupa fc Parmentie7 ( 2008[ ). who show the Milky- Way 
globular cluste r mass function to emerg e from a power-law embedded-cluster 
mass function. Parmentier et al. I ( 2008 ) expand on this by studying the effect 
that different assumptions on the physics of gas removal have on shaping the 
star-cluster mass function within about 50 Myr. 

The general ansatz that residual gas expulsion plays a dominant role in 
early cluster evolution may thus bear the solution to the long-standing prob- 
lem that the deduced initial cluster mass function needs to have this turnover, 
while the observed mass functions of young clusters are feature-less power-law 
distributions. 



The origin of population II stellar halos 



The above theory implies naturally that a major field-star component is gen- 
erated whenever a population of star clusters forms. About 12Gyr ago, the 
MW began its assembly by an initial burst of star formation throughout a 
volume spanning about lOkpc in radius. In this volume, the star formation 
rate must have re ached lOMpJyT such that star clu sters with masses up to 



10 M0 formed (jWeidner. Kr oupa fc Larsen 1120041 ). probably in a chaotic 



turbulent early interstellar medium. The vast majority of embedded clusters 
suffered infant weight loss or mortality, the surviving long-lived clusters evolv- 
ing to globular clusters. The so generated field population is the spheroidal 
population II halo, which has the same chemical properties as the surviving 
(globular) star clusters, apart from enrichment effects evident in the most 
massive clusters. All of these characterist i cs em erge naturally in the above 
model , as pointed out bv lKroupa fc Boilv I (|2002|). bvlParmentier fc Gilmore 



(|2007l ) and most recently bv lBaumgardt. Kroupa fc Parmentier I (j2008l ). 



1.3 Long term, or classical, cluster evolution 

The long-term evolution of star clusters that survive infant weight loss and 
the mass loss from evolving stars is characterised by three physical processes: 
the drive of the self-gravitating system towards energy equipartition, stellar 
evolution processes and the heating or forcing of the system through external 
tides. One emphasis of star-cluster work in this context is on testing stellar- 
evolution theory and on the interrelation of stellar astrophysics with stellar 



20 Pavel Kroupa 



dynamics given that the steUar-evolution and t he dynamical-evolution t ime- 
scales are comparable . The reader is directed to Mevlan fc Heggid ( 1997 ) and 
iHeggie fc H"ilt1 (|2003f) for further details. 



Tidal tails 



Tidal tails contain the stars evaporating from long-lived star clusters (the mov- 
ing group II above). The typical S-shaped morphology of tidal tails close to the 
cluster are easily understood: Stars that leave the cluster with a slightly higher 
galactocentric velocity than the cluster are on slightly outward directed galac- 
tic orbits and therefore fall behind the cluster as the angular velocity about 
the galactic centre decreases with distance. The outward directed trailing arm 
develops. Stars that leave the cluster with slower galactocentric velocities than 
the cluster fall towards the galaxy and overtake the cluster. 

Given that energy equipartition leads to a filtering in energy space of the 
stars that escape at a particular time, one expects a gradient in the stellar 
mass function progressing along a tidal tail towards the cluster such that the 
mass function becomes flatter, i.e. richer in more massive stars. This effect is 
difficu lt to detect, but for exam ple the long tidal tails found emanating from 
Pal 5 (jOdenkirchen et al. II2003I) may show evidence for this. 



As emphasised bv lOdenkirchen et al. I (|2003f ). tidal tails have another very 



interesting use: they probe the gravitational potential of the Milky Way if 
the differential motions along the tidal tail can be measured. They are thus 
important future tests of gravitational physics. 



Death and hierarchical multiple stellar systems 

Nothing lasts forever, and star clusters that survive initial re-virialisation after 
residual gas expulsion and mass loss from stellar evolution ultimately cease 
existing after evaporating all member stars leaving a binary or a long-lived 
highly hierarchical multiple system composed of near-equal mass components 
La Fuente Marcos I ll997ill998h . Note that these need not be stars. These 
cluster remnants are interesting, because they may ac count for most of the 



hiera rchical multiple stellar systems in the Galactic field (jGoodwin fc Kroupa 



20051) with the implication that they would not be a product of star formation 



but rather of star-cluster dynamics. 



1.4 What is a galaxy? 

Star clusters, dwarf-spheroidal (dSph) and dwarf-elliptical (dE) galaxies as 
well as galactic bulges and giant elliptical (E) galaxies are all stellar-dynamical 
systems that are supported by random stellar motions, i.e. they are pressure- 
supported. But why is one class of these pressure supported systems referred 
to as star clusters, while the others are galaxies? Is there some fundamental 
physical difference between these two classes of systems? 
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Considering the radius as a function of mass, it becomes apparent that 
systems with M ^10^ Mq do not show a mass-radius relation (MRR) and 
have r ~ 4pc. More massive objects, h o wever , show a well-defined MRR. In 
fact. iDabringhausen. Hilker fc Kroupal ( 20081 ) find that the massive compact 
objects (MCOs), which have 10^ < M/Mq < 10^ , lie on the MRR of giant E 
galaxies 10^'^ M0) down to normal E galaxies (lO^^^M©), as is evident in 
fig- 11 

/ i\ f \ 0.60±0.02 

E/pc=10-3.-(^) . (28) 




10"^ 10"'' 10° 10^ 10^ 10^ lO"^ 10^ 10^ 10^ 



dynamical mass [ 1 0^ solar units] 

Fig. 3. Mass-radius d ata in dependence of the dynamical mass of pressure-supported 
stellar systems (from [Pabringhausen. Hilker fc Kroupa II2OO8I I. MCOs are massive 
compact objects (also referred to as ultra compact dwarf galaxies). The solid and 
dashed lines refer to eq. 1281 while the dash-dotted line is a fit to dSph and dE 
galaxies. 



Noteworthy is that systems with M ?J 10® Mq also sport complex stel- 
lar populations, while less massive systems have single-age, sing le-metallicity 
populations. Remarkably. IPflamm-Altenburg fc Kroupa" ( 2008 ) show that a 
stellar system with iV ^ 10® Mq and a radius as observed for globular clusters 
can accrete gas from a co-moving warm inter-stellar medium and may re-start 
star formation. The median two-body relaxation time is longer than a Hubble 
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time for Af <; 3 x 10^ Mq, and only for these systems is there evidence for a 
shght increase in the dynamical mass-to-light ratio. Intriguingly, {M/L)y k, 2 
for M < 10^ Mq, while {M/L)v « 5 for M > 10^ Mq with a possible decrease 
for M > 10^ Mq (fig. S]). Finally, the average stellar density maximises at 
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□ dwarf spheroidal galaxies □ 

□ 




luminosity [ 1 solar units] 

Fig. 4. Dynamical M/L values in dependence of the V-band luminosity of pressure- 
supported stellar systems ffrom lPabringhausen. Hilker fc Kroupa 1 12008!). MCOs are 
massive compact objects (also referred to as ultra compact dwarf galaxies). 



M = 10^ Mq with about 3 x 10^ Mg/pc^ (jPabringhausen. Hilker fc Kroupa 
20081) . 



Thus, 

• the mass 10^ Mq appears to be special, 

• stellar populations become complex above this mass, 

• evidence for some dark matter only appears in systems that have a median 
two-body relaxation time longer than a Hubble time, 

• dSph galaxies are the only stellar-dynamical systems with 10 < [M / L)v < 
1000 and as such are total outliers. 

• 10^ Mq is a lower accretion limit for massive star clusters immersed in a 
warm inter-stellar medium. 

M « 10^ Mq therefore appears to be a characteristic mass scale such that 
less-massive objects show characteristics of star clusters being well-described 
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by Newtonian dynamics, while more massive objects show behaviour more 
typical of galaxies. Defining a galaxy as a stellar-dynamical object which has 
a median two-body relaxation time longer than a Hubble time, i.e. essentially 
a system with a smooth pote ntial, may be an objective and useful way to 
define a galaxy ((Kroupa lllQQSh . 

Why only smooth systems show evidence for dark matter remains at best 
a striking coincidence, at worst it may be symptomatic of a problem in un- 
derstanding dynamics in such systems. 
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2 Initial 6D conditions 

The previous section gave an outline of some of the issues at stake in the 
realm of pressure-supported stellar systems. In order to attack these and other 
problems, we need to know how to set-up such systems in the computer. 
Indeed, as much as analytical solutions may be preferred, the mathematical 
and physical complexities of dense stellar systems leave no alternatives other 
than to resort to full-scale numerical integration of the 6N coupled first-order 
differential equations that describe the motion of the system through 67V- 
dimcnsional phase-space. There are three related questions to ponder: Given 
a well-developed cluster, how is one to set it up in order to evolve it forward 
in time? How does a cluster form, and how does the formation process affect 
its later properties? How do we describe a realistic stellar population (IMF, 
binaries)? Each of these is dealt with in the following sections. 



2.1 6D structure of classical clusters 



Because the state of a star cluster is never known exactly, it is necessary to 
perform numerical experiments with conditions that are, statistically, consis- 
tent with the cluster snapshot. To ensure meaningful statistical results for 
systems with only a few stars, say N < 5000, many numerical renditions of 
the same object are thus necessary. For example, systems with N = 100 stars 
evolve erratically and numerical experiments are required to map out the 
range of possible states at a particular time: the range of half-mass radii at 
an age of 20 Myr in 1000 numerical experiments of a cluster initially with 
N — 100 stars and with an initial half-mass radius ro.s — 0.5 pc can be 
compared with an actually observed object for testing consistency with the 
ini tial condition s . Exc el lent recent examples o f this approach can be found 



Hurlev et al.l (12005 ): Portegies Zwart et al.l ( 2007 ). with a recent review 
available bv lHut et al.l (j2007l ). and two text books have been written dealing 



with co mputational and more general aspect s of the physics of dense stellar 
systems (|Aarseth Il20^ iHeggie fc Hut II2003I ). 

The six-dimensional structure of a pressure-supported stellar system at 
time t is conveniently described by the phase-space distribution function, 
f{r,v;t), where r,v are the phase-space variables, and 



dN = f{r,v;t)d^x(fv 



(29) 



is the number of stars in 6D phase-space volume element d'^xd^v. I n the 
case of a steady-state, the Jeans theorem (jBinnev fc Tremaind Il987l their 
chapter 4.4) allows us to express / in terms of the integrals of motion, i.e. the 
energy and angular momentum. The phase-space distribution function can 
then be written 

f = f{r,v) = f{e,,l), (30) 



where 



Dense Stellar Populations: Initial Conditions 



25 



e, = lv' + <Pir) (31) 
is the specific energy of a star, and 

/ = |r X ti| (32) 
is tlie specific orbital angular momentum of a star. The Poisson equation is 

\/^<P{r) = 4:nGpm{r)^4:TrG J mfSv, (33) 

or in spherical symmetry, 

where is the phase-space mass-density of all matter and is equal to to / 
for a system with equal-mass stars. Most pressure-supported systems have a 
near-spherical shape and so in most numerical work it is convenient to assume 
spherical symmetry. 

For convenience it is useful to introduce the relative potential!^, 

'P=-'P + $Q (35) 

and the relative energy 

£ = -ee+^o = ^-\v^, (36) 

where is a constant such that / > for £ > Q and / = for £ < 0. 
The Poisson equation becomes V^'I' = —AttG pm subject to the boundary 
condition ']/ — > as r — > oo. 

One important property of stellar systems is the anisotropy of their velocity 
distribution function. Defining the anisotropy parameter 

V 

where Vg , are the mean squared tangential and radial velocities at a par- 
ticular location r, respectively. It follows that systems with (3 = everywhere 
have an isotropic velocity distribution function. 

If / only depends on the energy then the mean squared radial and tan- 
gential velocities are, respectively. 



P Jail vol. 



dvr dve dV(j) (38) 



The following discussion is based on iBinnev fc Tremaind (|l987l ). 
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and 



1 



P Jail vol. 



dVr dvg dVcf, . 



(39) 



If the labels and r are exchanged in eq. [311 it can be seen that one arrives 
at eq. [351 Eq. [3S1 and [351 are thus identical, apart from the labeling. Thus if 
/ — f{£), P = and the velocity distribution function is isotropic. 

If / depends on the energy and the orbital angular momentum of the stars 
(|Z| = |r X v\) then the mean squared radial and tangential velocities are, 
respectively. 



— 1 



P Jail vol. 



1 



and 



1 



P J all vol. 



vef 



W--{vl+vl 



dvr dvg dv^ (40) 



dvrdvgdv^. (41) 



If the labels 9 and r are exchanged in eq. [411 it can be seen that this time 
one does not arrive at eq. [40l Thus if / = f{£, I) then (3^0 and the velocity 
distribution function is not isotropic. 

This serves to demonstrate an elementary but useful property of the phase- 
space distribution function. 

A very useful series of distribution functions can be arrived at from the 
following simple form: 



The mass density, 

p^{r) =A'kF 



FE' 



/2!f' 



f > 0, 
f < 0. 



(42) 



V dv, 



(43) 



where the upper integration bound is given by the escape condition, £ 
- (l/2)i;^ = 0. Substituting = 2if'cos26l for some 9 leads to 



Pmir) 



Cn 







If > 0, 
If < 0. 



(44) 



For Cn to be finite, n > 1/2, i.e. homogeneous (n — 0) systems are excluded. 

The Lane-Emden equation follows from the spherically symmetric Poisson 
equation after introducing dimensionless variables s = r/b,tp = S'/S'o, where 



b= (4 7rGif, 



n-l 



and tfo = •f(O), 



1 d 

ds 



ds 







^ <o. 



(45) 
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H. Lane and R. Emden worked with this equation in the context of self- 
gravitating polytropic gas spheres which have an equation of state 

P = Kpl,, (46) 

where is a constant and p the pressure. It can be shown that 7 = 1 + 

i.e. that the density distribution of a stellar polytrope of index n is the same 

as that of a polytropic gas sphere with index 7. 

The natural boundary conditions to be imposed on eq. US] are at s = 0, 

1. ?A = 1 since ^{Q) = tf^o, and 

2. dip/ds = because the gravitational force must vanish at the centre. 

Analytical solutions to the Lane-Emden equation are possible only for a 
few values of n, remembering that a homogeneous (n = 0) stellar density 
distribution has already been excluded as a viable solution of the general 
power-law phase-space distribution function. 



The Plummer model 

A particularly useful case is 



= , (47) 

l + is2 

It follows immediately that this is a solution of the Lane-Emden equation for 
n = 5, and it also satisfies the two boundary conditions above, and therewith 
constitutes a physically sensible potential. By integrating the Poisson equation 
it can be shown that the total mass of this distribution function is finite, 

M^^Vs^ob/G, (48) 

although the density distribution has no boundary. The distribution function 
is 



with the relative potential 



^= , ^ (50) 



and density law 



Pm,Q 
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with the above tota l mass. Thi s dens ity distribution is known as the Plummer 
model named after I Plummer I (jlQllI ) where he showed that the density dis- 



tribution resulting from this model provides a reasonable, and in particular 
very simple analytical description of globular clusters. The Plummer model 
is, in fact, a work-horse for many applications in stellar dynamics because 
many of its properties such as the projected velocity dispersion profile can be 
calculated analytically. Such formulae are useful for checking numerical codes 
used to set-up models of stellar systems. 



Properties of the Plummer model 



Some us eful analytical result s can be derived for the Plummer density law 
(see also lHeggie fc Hut Il2003l their p. 73, for another compilation). 

For the Plummer law of mass Mod the mass-density profile (eq. I5ip can 
be written 



Pm{r) 



3 AL, 



1 



(52) 



The central number density is thus 



3iV 

The mass within radius r follows from M{r) = An Pm{r') r'^ dr' , 



(53) 



Mir) = Afeci ■ 







1 + 







Thus, 

Tpi contains 35.4 per cent of the mass, 

2 Tpi contain 71.6 per cent, 

5 Tpi contain 94.3 per cent and 

lOrpi contain 98.5 per cent of the total mass. 

The half-mass radius contains 50 per cent of the mass. 



(23 - i)-2 « 1.305 r. 



pi- 



(54) 



(55) 



The projected surface mass density, Sm{R) = 2 Pm{r) dz, where R is 
the projected radial distance from the cluster centre and Z is the integration 
variable along the line-of-sight (r^ ^ + Z^), is 
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Assume there is no mass segregation so that the mass-to-light ratio, T = 
(M/L), measured in some photometric system is independent of radius. The 
integrated light within projected radius R is 



/(i?) = (l/T) / Sp{R')2TiR' dR\ 
Jo 



KR) 



Mod r. 



pi 



r 



1 



i?2 



(57) 



(58) 



'pi ^'pi_ 

Thus, Tpi is the half-light radius of the projected star cluster, I{rp\) = 
0.5/(cx)). 

In the above equations p{r) — p,n{r)/rn, N{r) = M{r)/m and i7„ = 
Sp/rri are, respectively, the stellar number density, the number of stars within 
radius r and the projected surface number density profile if there is no mass 
segregation within the cluster, the average stellar mass, m, therefore not being 
a function of radius. 

The velocity dispersion can be calculated at any radius from Jeans eq. 11201 
For an isotropic velocity distribution {a^ = ct^ = cr^), such as the Plummer 
model, the Jeans equation yields 



1 



p(r) 



■ GM{r') , 
p[r ) 2 ' 



(59) 



since d(l){r) jdr = GM (r) / , and the integration bounds have been chosen to 
make use of the vanishing Pm{i') as r — > oo. Note that the above equation is 
also valid if M(r) consists of more than one spherical component such as a 
distinct core plus an extended halo. Combining eqs [52l [54l and [59l leads to 



ccl 



2r 



pi 



1 



1 + — 



(60) 



where cr(r) is the three-dimensional velocity dispersion of the Plummer sphere 
at radius r, a'^{r) = X]fc=r e '''K'^) ^"^(t) = 3f7^j3(r) since isotropy is 
assumed. 



A star with mass m positioned at r and with speed v = (X]fc=i'^fc 



1/2 



can escape from the cluster if it has a total energy Cbind = Ckin + Cpot = 
0.5 TO -I- ■m(f){r) > such that v > Vosdr), implying for the escape speed at 
radius r, Vcsc{r) ~ \/2 \ (j){r)\. The potential at r is given by the mass within r 
plus the potential contributed by the surrounding matter which is calculated 
by integrating the contributions from each radial mass shell. 



G 



M{r 
r 

GAL, 



1 



+ G- p{r )4:nr'^dr 

I nr' 



rpi / [1 + (r/rpi) 



211/2 ■ 



(61) 
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SO that 

voscir) = [ r^l^ ) : (62) 

The circular speed, Vc, of a star moving on a circular orbit at a distance 
r from the cluster centre is obtained from centrifugal acceleration, v'^/r = 
d(t)[r)/dr = GM{r)/r^, 

.-(^)^,|^. 

In many but not all instances of interest the initial cluster model is chosen 
to be in the state of virial equilibrium. That is, the kinetic and potential 
energies of each star balance such that the whole cluster is stationary. The 
scalar virial theorem, 

2K + W = 0, (64) 
where K and W are the total kinetic and potential energy of the cluster^ 



1 f°° 

2 Jo 



— — for the Phmrmer sphere, (65) 

64 Tpi 

W—— / (f){r) p{r) Anr" dr ^ 
2 Jo 

_ ^ — cci ^-|_^^ Plummer sphere. (66) 
32 rpi 

The total, or binding, energy of the cluster, i^tot ^W + K^is 

Etot = -K^^W. (67) 



be defined as a^j = 2 K/Mcd so that 



The characteristic three-dimensional velocity dispersion of a cluster can 

.2 
cl 

o-ci = ^ , 68) 

32 Tpi 

G A/ncl /„„N 

(69) 



rav 



..'(^) (70, 

introducing the gravitational radius of the cluster, Tg^av = GM^^j/lVKj. For 
the Plummer sphere rg^av — (32/3 7r)rpi = 3.4rpi, and the structure factor 



Eq.3.251.4 on p. 295 in iGradshtevn fc RvzhikI (|l980l ) are useful to solve the inte- 
grals for the Plummer sphere 
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32 

« 0.88. (71) 

Defining the virial ratio, 

a cluster can initially be in three possible states: 

^ virial equilibrium, 

expanding, (73) 
collapsing. 

Note that if initially Q < 1/2 the value Q = 1/2 will be reached temporarily 
during collapse, after which Q increases further until the cluster set tles in 
virial equilibrium after this violent relaxation phase ( Binnev fc Tremainell987 . 
p. 271). 

The characteristic crossing time through the Plummer cluster, 

tcr ^ (74) 
CTlD,cl 




128\2 -1 3 



using the characteristic one-dimensional velocity dispersion, (Tid,c1 = ctci/V^- 
Observationally the core radius is that radius where the projected surface 
density falls to half its central value. For a real cluster it is much easier to 
determine than the other characteristic radii. For the Plummer sphere, 

-Rcorc = (V2 - l) ' rpi = 0.64 rpi, (76) 



from eq. [56] assuming the mass-to-light ratio, T, is independent of radius. For 
a King model 

«s".-r^^y, (77, 



V47rG pm{0) 

is the King radius. From eq. [501 "'^(O) = GMoci/(2rpi) and from eq. [5H 
p„(0) = 3Meci/(47rr3j) so that 
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The singular isothermal model 

Another useful set of distribution functions can be arrived at by consider- 
ing n = oo. The Lane-Emden equation is not well defined in this limit, 
but for a polytropic gas sphere eq. [46] implies 7 — > 1 as n — > 00. Thus 
p = K pm, which is the equation of state of an isothermal ideal gas with 
K ^ kB T/rUp, where /cb is Boltzmann's constant, T the temperature and 
mp the mass of a gas particle. From the equation of hydrostatic support, 
dp/dr = -pmiGM{r)/r'^), where M(r) is the mass within r, the following 
equation can be derived 

dr [ dr J ^ kBT ^™ ^'^> 

For a distribution function (our ansatz) 

/™(f) = ^^e*, (80) 
(27r cr2)^ 

where is a new quantity related to a velocity dispersion and remembering 
8 = 9 — i'^/2, one obtains from pm = / fm{£) 4 7ru^ 

^^r)=ln(^^] a\ (81) 

From the Poisson equation it then follows that 

a = const = ^e2- (32) 
mp 

for consistency with eq. [751 

Therefore, the structure of an isothermal, self-gravitating ideal sphere of 
gas is identical to the structure of a coUisionless system of stars whose phase- 
space mass-density distribution function is given by eq. [8OI Note that f{£) is 
non-zero at all £ (cf to King's models below) . 

The number- distribution function of velocities is F{v) = J^^^ f{£)d^x, 
i.e. ^ 

F{v)=Fae-f^. (83) 

This is the Maxwell-Boltzmann distribution which results from the kinetic 
theory of atoms in a gas at temperature T that are allowed to bounce off 
each other elastically. This exact correspondence between a stellar-dynamical 
system and a gaseous polytrope holds only for an isothermal case (n = 00). 

The total number of stars in the system is iVtot = -/Vtot /q°° F{v) iTTv"^ dv 
and the number of stars in the speed interval v to v + dv is 



dN = F{v)'iTrv^ dv = Ntot -e ^Airv^dv, (84) 

(2 7rcr2)2 
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which is the Maxwell-Boltzmann distribution of speeds. The mean-squared 
speed of stars at any point in the isothermal sphere is 

— 47r v"^ F{v)v'^ dv 



47r F{v)v'^dv = 3fT2' 

and the one-dimensional velocity dispersion is ctid = ctq = ct, where a — 
r,e,(t),x,y,z, .... 

To obtain the radial mass-density of this model, the ansatz p,n — C r^'' 
together with Poisson's equation (eq. [79|) implies 

Pm{r)^^^, (85) 
i.e. a singular isothermal sphere. 

The isothermal model 

The above model has a singularity at the origin, which is unphysical. In order 
to remove this problem, it is possible to force the central density to be finite. To 
this end new dimensionless variables are introduced, p'm = Pm/ Pmfi, f = rjr^. 
p~m is the finite central density, while ro = P^orc is the King radius (eq. [77|) 
such that the projected density falls to 0.5013 (i.e. about half) its central 
value. To is also sometimes called the core radius (but see further below for 
King models on p. [35|). Poisson's equation (eq. [79]) therewith becomes 

YrK ^)=-^<'-' ■ 

This differential equation must be solved numerically for Pm(f) subject to the 
boundary conditions (as before), 

p™(r-0) = l, ^1,^0-0. (87) 
dr 

The solution is the isothermal sphere. 

By imposing physical reality (central non-singularity) onto our mathemat- 
ical ansatz we end-up with a density profile which cannot be arrived at analyt- 
ically but only numerically. The isothermal density sphere must be tabulated 
in the computer with entries such as 

r/ro, logio ( — ) logio ^ 



PoJ \ro po 

where S is the pr ojected density (for example see table 4-1 and fig. 4-7 in 



Binnev fc Tremain c 1987). The circular velocity, Vc{r) — GM(r)/r of the 



isothermal sphere obtains by integrating Poisson's equation (eq. [79)1 from r = 
to r = r' such that r"^ {dlnp„J dr) = -{G/a"^) M{r) and 
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2 2 d\np„, {r) 

Numerical solution of differential eg . [551 shows that Vc — > \/2(t (constant) 
for large r. 

The isothermal sphere is a useful model for describing elliptical galaxies 
within a few core radii and disk galaxies because of the constant rotation 
curve. However, combining the two equations for above, one finds that 
M(r) « (2a^/G)r for large r, i.e. the isothermal sphere has an infinite mass 
as it is not bounded. 



The lowered isothermal or King model 



We have thus seen that the class of models with n = oo contains as the simplest 
case the singular isothermal sphere. By forcing the central density to be finite 
we are led to the isothermal sphere, which however, has an infinite mass. The 
final model considered here within this class is the lowered isothermal model, 
or the King mode@, which forces not only a finite central density but also 
a cutoff in radius. These have a distribution function similar to that of the 
isothermal model, except for a cutoff in energy, 



(2 TTCr^) 2 



- Ij : £>0, 
: £<0. 



(90) 



The density distribution becomes 



PrnA 



erf 



4if 



1 



2!f 
3^ 



(91) 



integrating only to £ = as before. The Poisson eq. [72] becomes 



d / 2 dlnp'j, 
df \ df 



erf 



4^ [ 



(92) 

Again, this differential equation must be solved numerically for !f'(r) subject 
to the boundary conditions. 



(93) 



Note that [Kind (| 19621 ) suggested three-parameter (mass, core radius and cut- 
off/tidal radius) empirical projected (2D) density laws that fit globular clusters 
very well. These do not have information on the velocity structure of the clusters. 
The King (non-analytical) 6D models that are solu tions o f the Jeans equation 
feg. 11201 below) and discussed here, are published by [Kind l|l966l ). 
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The density vanishes at r = 7-tid (the tidal radius), where 'F{r — rtid) = 
also. A King model is thus limited in mass, has a finite central density but the 
parameter a is not the velocity dispersion but is rather related to the depth 
of the potential via the concentration parameter 

. £<5>. (94) 

The concentration is defined as 

c^log,o(^^). (95) 

For globular clusters, 3 < Wo < 9, 0.75 < c < 1.75, and the relation between 
Wo and c is plotted in fig. [5] Note also that the true core radius defined as 



O 1.5 




Wo 

F ig. 5. The KinR-conc entration parameter Wo as a function of c (cf. with fig. 4-10 
in lBinnev fc TremainellT987i ). This figure has been produced by Andreas Kiipper. 



S{Rc) = (1/2) ^(0), where S{R) is the projected density profile and R is the 
projected radius, is unequal in general to the King radius, ro fea. [77]). 

Finally, it should be emphasised that it is not physical to use an arbitrary 
rtid^ the tidal radius must always match the value dictated by the cluster mass 
and the host galaxy (e.g. eq. [3]). 

2.2 Comparison: Plummer vs King models 

The above discussion has served to show how various popular models can 
be followed through from a power-law distribution function (eq. 1421) with 
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different indices n. The Plummer model (p. [77)1 and the King model (p. [M]) 
are particularly useful for describing star clusters. The Plummer model is 
determined by two parameters, the mass, M, and the scale radius, tq.s ~ 
1.305 Tpi. The King model requires three parameters, M, a scale radius, rp.s, 
and a concentration parameter. Wo or c. Which sub-set of parameters yield 
models that are similar in terms of the overall density profile? 

To answer this, the mass is set to be constant. King models with different 
Wo and ro.5 are computed and Plummer models are sought that minimise the 
reduced chi-square value between the two density profiles. Fig. [S] shows two 
examples of best-matching density profiles, and fig. [7] uncovers the family of 
Plummer profiles that best match King models with different concentration. 
Note that a good match between the two is only obtained for intermediate 
concentration King models (2.5 ^ Wo ~ 7.5). 



2.3 Discretisation 

To set up a computer model of a stellar system with N particles (e.g. stars) the 
distribution functions need to be sampled N times. The relevant distribution 
functions are the phase-space distribution function, the stellar initial mass 
function and the three distribution functions governing the properties of binary 
stars (periods, mass-ratios, eccentricities). 

Assume the distribution function depends on the variable Cmin < C < Cmax 
(e.g. stell ar mass, m). Ther e are various ways of sampling from a distribution 
function ( Press et al. 1992f) . but the most efficient way is to use a generating 



function if one exists. Consider the probability, X{(), of encountering a value 
for the variable in the range Cmin to C, 

XiO^I^^ piOdC, (96) 

with X{(^ynin) = < X{Q < X{(^ynax) = 1, and p{C) is the distribution 
function normalised such that the latter equal sign holds {X = 1). p{C,) is the 
probability density. The inverse of eq. 1961 C(^)j is the generating function. It 
is a one-to-one map of the uniform distribution X G [0, 1] to C G [Cmin, Cmax]. 
If an analytical inverse does not exist, then it can be found numerically in 
a straight-forward manner for example by constructing a table of X, C and 
then interpolating this table to obtain a C, for a given X. 

Example: The power-law stellar mass function 

As an example, consider the distribution function 

Tfl 

^{m) = km-'^, a = 2.35; 0.5 < — < 150. (97) 

The corresponding probability density is p{m) = kprn""', and J^^^ p{m) dm — 
1 fcp 0.53. Thus 
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r/r^ 

Fig. 6. Comparison of a King model (solid curve) with a Plummer model (dashed 
curve). Both have the same mass, and that Plummer model is sought which min- 
imizes the unweighted reduced chi-squared between the two models. The upper 
panel shows a high-concentration King model with c = 2.55, Wo ~ 11, and the 
best-fit Plummer model has rj's""™'^'^ = 0.366 r^™" {rh = '"o.s), as stated in the 
panel. The lower panel compares the two best matching models for the case of an 
intermediate-concentration King model. This figure has been produced by Andreas 
Kiipper. 
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Fig. 7. The ratio r^Sr^"""' /r^'^'^ {rn = ro.s) for tlie best-matcliing Plummer and 
King models (fig. |6]) are plotted as a function of the King concentration parameter 
Wo- The uncertainties are unweighted reduced chi-squared values between the two 
density profiles. It is evident that there are no well-matching Plummer models for 
low- (c < 2.5) and high-concentration (c > 7.5) King models. This figure has been 
produced by Andreas Kiipper. 



X{m) 



p{m) dm 



150 



0.5^ 



1 



/0.5 

and the generating function for stellar masses becomes 



m(X) 



X- 



1 



a 



0.5 



l-a 



(99) 



It is easy to programme this into an algorithm: obtain a random variate X 
from a random number generator and use the above generating function to 
get a corresponding mass, m. Repeat N times. 



Generating a Plummer model 

Perhaps the most useful and simplest model of a bound stellar system is 
the Plummer model (p. Wf\ . It is worth introducing the discretisation of this 
model in some detail, because analytical formulae go a long way which is 
important for testing codes. A condensed form of this material is available in 
Aarseth et all (|l974l ). 

The mass within radius r is (7'pi = b here) 



M{r)= / p„,(/)4^r'2dr' = Mci 
Jo 



1 + {r/r^x)' 



(100) 
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A number uniformly distributed between zero and one can then be defined, 

where ( = r/r-p\ and Xi{r = oo) = 1. This function can be inverted to 
yield the generating function for particle distances distributed according to a 
Plummer density law, 

_ 1^ 

c(Xi) = (xr* - 1) ^ (102) 

The coordinates of the particles, x, y, z, = (C^'pi)'^ — + + z^, can be 
obtained as follows: For a given particle we already have r. For all possible 
X and y, z has a uniform distribution, p{z) = const = 1/(2 r) over the range 
— r < z < +r. Thus, for a second random variate between zero and one, 

X2{z) = y p{z') dz' =l-{z + r), (103) 

with X2{+r) = 1. The generating function for z becomes 

z{X2) = 2rX2-r. (104) 

Having obtained r and z, x and y can be arrived at as follows, noting the 
equation for a circle, r"^ — z^ = a;^ + y^: Choose a random angle 6 which is 
uniformly distributed over the range < ^ < 27r. Thus p{6) = 1/(2 tt) and 
the third random variate becomes 

X3{e)= f ^de' = ^. (105) 

The corresponding generating function is 

^(X3)=2 7rX3. (106) 

Finally, 

x = {r^ - z^) ^ cosO; y = {r^ - z^) ^ sinO. (107) 

The velocity for each particle cannot be obtained as simply as the positions. 
In order for the initial stellar system to be in virial equilibrium, the potential 
and kinetic energy need to balance according to the scalar virial theorem. 
This is ensured by forcing the velocity distribution function to be that of the 
Plummer model, 

/„(ee) = I • (108) 

[ : Ce > 0, 



where 
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ee(r,w) =^>(r) + (l/2)w2 (109) 
is the specific energy per star, and 



is the potential. Now, the Plummer distribution function can be expressed in 
terms of r and w, 

f{r,v)^h(^-^r)-]^v^y , (111) 

for a normalisation constant fo and dropping the mass sub-script because we 
assume the positions and velocities do not depend on particle mass. With the 
escape speed at distance r from the Plummer centre, Vescir) = a/— 2 <l>{r) = 
ii/C, it follows that 

f{r,v)^fo Qi-csc) {i-ey ■ (112) 

The number of particles with speeds in the interval u, u + is 

dN ^ f{r,v)'iTTv^ dv = g{y)dv. (113) 

Thus g 

5W = 16 7rL QT;esc(r)) (l - C'M) ^ C'(r), (114) 

that is, 

9{C,) = 9oC{r){l-e{r))\ (115) 
for a normalisation constant go determined by demanding that 



X4(C-l) = l= g{C)dC 
Jo 



(116) 



for a fourth random number deviate ^4(C) = Jq giC) dC ■ It follows that 

X4(C) = i (5C^-3C^). (117) 

This cannot be inverted to obtain an analytical generation function for C = 
({X4). Therefore, numerical methods need to be used to solve eg. 11171 For 
example, one way to obtain for a given random variate X4 is to find the 
root of the equati on = (1/2) (5 C^ - 3 C^) - ^4, or one can use the Neumann 
rejection method (jPress et al.lll992l ). 

The following procedure can be implemented to calculate the velocity vec- 
tor of a particle for which r and C are already known from above: Compute 
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^^csc(^) so that V = C'^^csc- Each speed v is then spht into its components 
VxjVy^v^ assuming velocity isotropy using the same algorithm as above for 
X J y , z'. 

v,{X5)^{2X5-l)v; 0{Xe)^27rXe; (118) 



Vx — \/ — cos9; Vy — y'v^ — v^sinO. (119) 

Note that a rotating Plummer model can be generated by simply switching 
the signs of Vx and Vy such that aU particles have the same direction of motion 
in the x — y plane. 

As an aside, an efficient numerical method to set-up triaxial spheroids with 
or without an embedded rotating disk is described bv lBoilv et all (|200l[ ). 



Generating an arbitrary spherical, non-rotating model 

In most cases an analytical density distribution is not known (e.g. the 
King models above). Such numerical models can nevertheless be discretised 
straightforwardly as follows. Assume that the density distribution, p{r), is 
known. Compute M(r) and M^. Define X{r) = M{r)/Mci, as above. We thus 
have a numerical grid of numbers r, M(r), X(r). For a given random deviate 
X E [0, 1], interpolate r from this grid. Compute above. 

If the distribution function of speeds is too complex to yield an analytical 
generating function X{() for the speeds C, then one can resort to the following 
procedure: One of the Jeans equations for a spherical system is 

^ (p(r)a.(r)2) + ^ [2aUr) - (a,(r)^ + ..(r)^)] = -p(r) (120) 

For velocity isotropy, = erg = this reduces to 

^ = -p?. (121) 
dr dr 

Integrating this by making use of p — > as r — > oo, and remembering that 
d$/dr = -GM/r^, 

For each particle at distance r a one-dimensional velocity dispersion, ar (r) , is 
thus obtained. Choosing randomly from a Gaussian distribution with disper- 
sion (Ti,i = r,9,(j),x,y, z then gives the velocity components (e.g. Vx,Vy,Vz) 
for this particle. 
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Rotating models 

Star clusters are probably born with some rotation because the pre-cluster 
cloud core is likely to have contracted from a cloud region with differential 
motions that do not cancel out. How large this initial angular momentum 
content of an embedded cluster is remains uncertain, as the dominant motions 
are random chaotic ones owing to the turbulent velocity field of the gas. Once 
the star formation process is quenched as a result of gas blow-out (section [TTT]) 
the cluster expands which must imply substantial reduction in the rotational 
velocity. A case in point i s uj Cen, which has b een found to rotate with a peak 
velocity of about 7km/s ( Pancino et al.|[2007 and references therein). 



Setting-up rotating cluster models is easily done, e.g. by increasing the tan- 
gential velocities of stars by a certain factor. A systematic study of relaxation- 
driven angular momentum re-distribution within star clusters has become 
available through the work of the group of Rainer Spurzem and Hyung-Mok 



Lee, and the interested reader is directed to that body of work ijKim et al 



20081 and references therein) . One important outcome of this work is that core 



collapse is substantially sped- up in rotating models. The primary reason for 
this is that increased rotational support reduces the role of support through 
random velocities for the same cluster dimension. Thus, the relative stellar 
velocities decrease and the stars exchange momentum and energy more effi- 
ciently, enhancing two-body relaxation and therewith the dive towards energy 
equipartition. 



2.4 Cluster birth and young clusters 

Some astro-physical issues related to the initial conditions of star clusters have 
been raised in section 11.11 In order to address most of these issues numerical 
experiments are required. The very initial phase, the first 0.5 Myr, can only be 
treated through gas-dynamical computations that, however, lack the numer- 
ical resolution for the high-precision stellar-dynamical integrations which are 
the essence of collisional dynamics during the gas-free phase of a cluster's life. 
This gas-free stage sets-in with the blow-out of residual gas at an age of about 
0.5-1.5 Myr. This time is dominated by the physics of stellar feedback and ra- 
diation transport in the residual gas as well as energy and momentum transfer 
to it through stellar outflows. The gas-dynamical computations cannot treat 
all the physical details of the processes acting during this critical time, which 
also includes early stellar-dynamical processes such as mass segregation and 
binary-binary encounters. 

One successful procedure to investigate the dominant macroscopic physical 
processes of these stellar-dynamical processes, gas blow-out and the ensuing 
cluster expansion through to the long-term evolution of the remnant cluster, 
is to approximate the residual gas component as a time-varying potential in 
which the young stellar population is trapped. The pioneering work using 
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this approach has bee n performed by Lada et al. ( 1984 ). whereby the ear- 
lier numerical w ork by Tutukov _ (ll978) on open clusters and later iV— body 
computations by lGoodwinl ( 1997al b. 1998) on globular clusters must also be 
mentioned in this context. 

The physical key quantities that govern the emergence of embedded clus- 
ters from their clouds and their subseque nt appearance are (section [TTT] and 
Baumgardt. Kroupa fc Parmentier 1120081 ): 



• sub-structuring; 

• initial mass segregation; 

• the dynamical state at feedback termination (dynamical equilibrium?, col- 
lapsing? or expanding?); 

• the star-formation efficiency, e; 

• the ratio of the gas-expulsion time-scale to the stellar crossing time through 
the embedded cluster, Tgas /Across; 

• the ratio of the embedded-cluster half-mass radius to its tidal radius, rh/rt. 

It becomes rather apparent that the physical processes governing the 
emergence of star clusters from their natal clouds is terribly messy, and the 
research-field is clearly observationally driven. Observations have shown that 
star clusters suffer substantial infant weight loss and probably about 90 per 
cent of all clusters disperse altogether (infant mortality) . This result is consis- 
tent with the observational insight that clusters form in a compact configura- 
tion with a low star formation efficiency (0.2 J; 0.4) and that residual-gas 
blow-out occurs on a time-s cale comparabl e or even faster than an embedded- 
cluster crossing time-scale ( Kroupa Il2005l) . Theoretical work can give a rea- 
sonable description of these empirical findings by compactifying some of the 
above parameters, such as working with an effective star- formation efficiency 
as being a measure of the amount of gas removed for a cluster of a given 
stellar mass assuming this cluster was in dynamical equilibrium at feedback 
termination, and that the gas and stars were distributed according to the 
same radial density function with the same scaling radius. 



Embedded clusters: One way of parametrising an embedded cluster is to 
set-up a Plummer model in which the stellar positions follow a density law 
with the parameters Mod and rpi, and the residual gas is a time- varying 
Plummer potential initially with the parameters Aigas and r^\ (i.e. same 
radial density law). The effective star-formation efhciency is then given by 
eq. [H Stellar velocities must then be calculated from a Plummer law with 
total mass Mod + Mgas following the recipes of section | 2.3I The gas can be 



remov ed b y evolving Mgas or r-pi. For example , iKroupa. Aarseth fc Hurlev 



(|200ll) and iBaumgardt. Kroupa fc Parmentier I (j2008[ ) assumed the gas- mass 
decreases on an exponential time-scale after an embedded phase lasting about 
0.5 Myr during which the cluster i s allowed to evolve in dynamical equi- 
librium. Bastian fc Goodwin ( 2006f ). as another example, do not include a 
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gas potential but take the initial velocities of stars to be '^j \f(. times larger, 
"ycmbcddcd = gas, to modcl the effect of instanta neous ga s remo val. 

Many variations of these assumptions are possible, and lAdamsl (|200nf ). for 



example, investigated the fraction of stars left in a cluster remnant if the 
radial scale length of the gas is different to that of the stars, i.e. for a radially- 
dependent star- formation efficiency, e(r). 



Sub-c lustering: Initial sub-clustering has been barely studied: IScallv fc Clarke 



( 2OO2I ) considered the degr ee of sub-structuri i ig of t he ONC allowed by its cur- 



rent morphology, while Fc Uhauer fc Kroupa ( 20051 ) computed the evolution of 



massive star-cluster complexes assuming each me mber cluster in the complex 



under goes its own individual gas-expulsion process. iMcMillan. Versperini fc Portegies Zwart 



(|2007l ) showed that initially mass-segregated sub-clusters retain mass-segregation 
upon merging, this being an interesting mechanism for speeding-up dynamical 
mass segregation as it occurs faster in smaller- iV systems which have a shorter 
relaxation time. 

The simplest initial conditions for such numerical experiments are to set-up 
the star-cluster complex (or proto-ONC-type cluster for example) as a Plum- 
mer model, where each particle is a smaller sub-cluster. Each sub-cluster is also 
a Plummer model, embedded in a gas potential given as a Plummer model. 
The gas-expulsion process from each sub-cluster can be treated as above. 

Mass segregation and gas blow-out: The problem of how initially mass- 
segregated clusters react to gas blow-out has not been studied at all in the 
past owing partially to the lack of convenient algorithms for setting-up mass- 
segregated clusters that are in dynamical equilibrium and which do not go 
into core-collapse as soon as the TV— body integration begins. An interesting 
aspect here is that gas blow-out will unbind mostly the low-mass stars, while 
the massive stars are retained. These, however, evolve rapidly such that the 
mass- lost from the remnant cluster owing to the evolution of the massive stars 
can become destructive, enhancing infant mortality. 

Ladislav Subr has developed a numerically efficient method for setting- 
up initially mass-segregated clusters close to core-collapse based on a novel 
concept using the potentials of s ub-sets of stars ordered by their mass 
( Subr. Kroupa fc Baumgardt 2008 ^. An alternative algorithm based on or- 



dering the stars by increasing mass and increasing total energy leading to 
total mass segregation, but also to a model that is not in core collapse 
and which therefore evolves towa r ds co re collapse, has been developed by 



Baumgardt. Kroupa fc de Marchi I (|2008[ ) 



An application concerning the effect on the observed stellar mass function 
in globular clusters shows that gas-expulsion leads to bottom-light stellar mass 



The C-language software package plumix can be down-loaded at 
jhttp : //www, astro .vini-bonn. de/~webaiub/english/downloads .phpl . 
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functi ons in clusters with a low concentra tion, consistent with observational 
data ( Marks. Kroupa fc Baumgardt 2008f ). 
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3 The stellar IMF 

The stellar initial mass function (IMF), ^(m) dm, where m is the stellar mass, 
is the parent distribution function of the masses of stars formed in one event. 
Here, the number of stars in the mass interval m, m + dm is 



dN = ^(m) dm. 



(123) 



The IMF is, strictly speaking, an abstract theoretical construct because any 
observed system of N stars merely constitutes a particular represento^jon of 
this universa l distribution function, if such a function exists ( Elmegreen II1997I: 
[m aiz Apellaniz fc Ubeda 1 12005) . The probable existence of a unique ^ (m) can 
be inferred f rom observatio ns of an ensemble of systems each consisting of 
N stars (e.g. ^Ma ssev 20031 ). If, after corrections for (a) stellar evolution, (b) 
unknown multiple stellar systems, and (c) stellar-dynamical biases, the indi- 
vidual distributions of stellar masses are similar within the expected statistical 
scatter, then we (the community) deduce that the hypothesis that the stellar 
mass distributions are not the same can be excluded. That is, we make the 
case for a universal, standard or canonical stellar IMF within the physical 
conditions probed by the relevant physical parameters (metallicity, density, 
mass) of the populations at hand. 

Related overviews of the IMF can be found in lKrour)al(|2002al ): IChabr ler 



(|2003f) : iBonnell et all (|2007t ): iKroupal (|2007a[) . and a review with an em 



phasis on the met a l-rich problem is available in iKroupal (|2007bh . while 



Zinnecker fc Yorke I (|2007l ) p rovide an in-depth review of the formation and 



distribution of massive stars. lElmegreenI (|2007l ) discusses the possibility that 
star-formation occurs in different modes with different IMFs. 



3.1 The canonical or standard form of the stellar IMF 

The canonical stellar IMF is a two-part-power law fea. ll28[) . the only structure 
with confidence found so far being the change of index from the Salpeter/Massey 
value to a smaller one near 0.5 Mg0: 



^(m) cx m 



1,2 



(124) 



ai = 1.3±0.3, , 0.08 < to/Mq< 0.5, 
a2 = 2.3±0.5, , 0.5< m/Mo<m^ax, 

where mmax < 'Tiniax* ~ 150 Mq follows from fig. [1] Brown dwarfs have 
been found to form a separate population with aa w 0.3 ± 0.5 (eq. 11291 
Thies fc Kro"upall2007t ). 



^° The uncer tainties in Oj are estimated from the alpha-plot (section |3.2|I . as shown 
in fis. 5 in iKroupa I l|2002bD . to be about 95 per cent confidence limits 
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It has been corrected for bias throug h unresolved multi ple stellar sys- 
tems in the low- mass (m < 1 Mq) regime ([Kroupa et al.lll991l) using a multi- 
dimensi onal optimisation t echnique. The general outline of this technique is as 
follows (jKroupa. Tout fc G ilmorc 1993): first the correct form of the stellar- 
mass-luminosity relation is extracted using observed stellar binaries and the- 
oretical constraints on the location, amplitude and shape of the minimum 
of its derivative, dm/dMy, near m ~ 0.3 M©, My ft! 12, M/ « 9 in combi- 
nation with the observed shape of the nearby and deep Galactic-field stellar 
luminosity function (LF) 



( dm 



^(m). 



(125) 



where dN = i?'(My) dMy is the number of stars in the magnitude inteval 
My to My + dMv- Having established the semi-empirical mass-luminosity 
relation of stars, which is an excellent fit to the most recent observational 



constraints bv lDelfosse et all (|2000l ). a model of the Galactic field is then cal- 
culated assuming a parametrised form for the MF and different values for the 
scale- height of the Galactic-disk, and different binary fractions in it. Measure- 
ment uncertainties and age and metallicity spreads must also be considered in 
the theoretical stellar population. Optimisation in this multi-parameter space 
(MF parameters, scale-height and binary population) against observational 
data leads to the canonical stellar MF for m < 1 Mq . 

One important result from this work is the finding that the LF of main 
sequence stars has a universal sharp peak near My ~ 12, Mj » 9. It results 
from changes in the internal constitution of stars that drive a non-linearity in 
the stellar mass-luminosity relation. 

A consistency-check is then performed as follows: the above MF is used in 
creating young populations of binary systems (section 14.21) that are born in 
modest star clusters consisting of a few hundred stars. Their dissolution into 
the Galactic field is computed with an A^-body code, and the resulting theo- 
retical field is compared to the observed LFs (fig. [8]). Further confirmation of 
th e form of the canon ical IMF comes from indep endent sources, most notably 
bv lReid et al"! (|2002h and al so IChabrier I (120031) . 

In the high-mass regime, iMassev I (|2003l ) reports the same slope or index 
as = 2.3 ± 0.1 for toJJIOMq in many OB associations and star clusters 
in the Milky Way (MW), the Large- and Small-Magellanic clouds (LMC, 
SMC, respectively). It is therefore suggested to refer to a2 = = 2.3 as the 
Salpeter/Massey slope or index, given the pioneering work of Salpeterl ( 1955h 
who derived this value for stars with masses 0.4 — 10 Mq. 

Multiplicity corrections await to be published once we learn more about 
how the components are distributed in massi ve stars fcf. lPreibisch et ab 1l999l : 
Zinneckedliool . lWeidner fc Kroupal (|2008l ) are in the process of performing 
a very detailed study of the influence of unresolved binary and higher-order 
multiple stars on determinations of the high-mass IMF. 
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Single stars 




Fig. 8. The Galactic-field population as resulting from disrupted star clusters: Uni- 
fication of both, the nearby (solid blue histogramme) and deep (filled red circles) 
LPs with one parent MP (eq. I124p . The theoretical nearby LP (the solid line) is 
simply the LP of all individual stars, while the dashed line is a theoretical LP with 
a mixture of about 50 per cent unresolved binaries and single stars stemming from 
a clustered star-formation mode. According to this model all stars are formed as bi- 
naries in modest clusters which disperse to the field, and the resulting Galactic-field 
population has a binary fraction and a mass-ratio distribution as observed. After 
iKroupal (Il995al lb[). Note the peak in both theoretical LPs - it stems from the ex- 
tremum in th e derivative of t he stellar mass-luminosity relation in the mass range 
0.2 - 0.4 Mq (lKroupall2002bl ). 



Contrary to the Salpeter/Massey index (a = 2.3) IScalo ( 19861 ) found 
^MWdisk ~ 2.7 (m^ 1 Mo) from a very thorough analysi s of OB star counts 
in the MW disk. Similarly, the star-count ana l ysis of Reid et al~l ( 2002 ) 
leads to 2.5 < aMWHi.v < 2. 8. and iTinslevI (Il980l) . iKennicuttI (Il983f) (his ex - 
tended Miller-Scalo IMF), iPortinari et all (|2004f) and iRomano et al.l (|2005l ) 
find 2.5 QiMWdisk ~ 2.7. That aMWdisk > 0^2 follows naturally is shown in 
section 13.41 

Below the hydrogen-burning limit (see also section [3?3)) there is substantial 
evidence that the IMF flatteiis furt her to ao « 0.3 ± 0.5 (jMartm et al.ll2000l : 
Chabrier|[2003l : lMoraux et al.ll2004l ). Therefore, the canonical IMF most Hkely 



has a peak at 0.08 Mq. Brown dwarfs, however, comprise only a few per cent 
of the mass of a population and are therefore dynamically irrelevant (tabled]). 
The logarithmic form of the canonical IMF, 
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CL(m) ^ In(lO) m ^(m), (126) 

which gives the number of stars in logiow-intervals, also has a peak near 
0.08 Mq. However, the system IMF (of stellar single and multiple systems 
combined to system masses) has a maximum in the mass range 0.4 — 0.6 



(jKroupa et al.l 120031 ) 



The above canonical or standard form has been derived from detailed 
considerations of Galactic-field star-counts thereby representing an average 
IMF: for low-mass stars it is a mixture of stellar populations spanning a large 
range of ages (0 — lOGyr) and metallicities ([Fe/H] ^ — 1). For the massive 
stars it constitutes a mixture of different metallicities ([Fe/H] ^ — 1.5) and 
star-forming conditions (OB associations to very dense star-burst clusters: 
R136 in the LMC). Therefore it can be taken as a canonical form, and the 
aim is to test the 



IMF Universality Hypothesis: the canonical IMF constitutes the parent 
distribution of all stellar populations. 



Negatio n of this hypo thesis would imply a variable IMF. Note that the 



work of iMassev I (|2003f ) has already established the IMF to be invariable for 
m^lOM© and for densities 10^ stars/pc^ and metallicity 0.002. 

Finally, table [2] compiles some numbers that are useful for simple insights 
into stellar populations. 



3.2 Universality of the IMF: resolved populations 

The strongest test of the IMF Universality Hypothesis is obtained by 
studying populations that can be resolved into individual stars. Since one also 
seeks co-eval populations with stars at the same distance and with the same 
metallicity to minimise uncertainties, star clusters and stellar associations 
would seem to be the test objects of choice. But before contemplating such 
work some lessons from stellar dynamics are useful: 



Star clusters and associations 



To access a pristine population one would consider observing star-clusters 
that are younger than a few Myr. However, such objects carry rather mas- 
sive di sadvantages: the pre-main sequence stellar evolutio n tracks are unre- 
liable (jBaraffe et al.ll2002t IWuchterl fc Tscharnuter]|2003l ) such that the de- 
rived masses are uncertain by at least a factor of about two. Remaining gas 
and dust lead to patchy obscuration. Very young clusters evolve rapidly: the 
dynamical crossing time is given by eq. |4] where the cluster radii are typically 
ro.5 < Ipc and for pre-cluster cloud-core masses Afgas+stars > 10'^ M© the 
velocity dispersion CTcI > 2km/s such that ta < 1 Myr. 

The inner regions of populous clusters have icr ~ 0.1 Myr, and thus sig- 
nificant mixing and relaxation occurs there by the time the residual gas is 
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mass range 




riN 






Vm 








[Mq] 


[per cent] 


[per cent] 




[Mq/pc-^] 


[Mq/pc^] 






as 






as 




as 


as 




2.3 


2.7 


4.5 


2.3 


2.7 


4.5 


4.5 


4.5 


n n 1 n no 


37.2 


37.7 


38.6 


4.1 


5.4 


7.4 


3.2 X 10~^ 


1.60 


0.08-0.5 


47.8 


48.5 


49.7 


26.6 


35.2 


48.2 


2.1 X 10"-^ 


10.5 


0.5-1 


8.9 


9.1 


9.3 


16.1 


21.3 


29.2 


1.3 X 10"^ 


6.4 


1-8 


5.7 


4.6 


2.4 


32.4 


30.3 


15.1 


6.5 X 10"^ 


1.2 


8 - 120 


0.4 


0.1 


0.0 


20.8 


7.8 


0.1 


3.6 X 10-^ 


6.5 X 10"^ 


m/MQ = 


0.38 


0.29 


0.22 




pll, = 0.043 


^tot = 19.6 






Q3 = 


= 2.3 


as - 


= 2.7 




ziMci/Mci 








Mci 


N,i 


Mci 


mto 


[per cent] 




[Mq] 








[Me] 




as = 2.3 


as = 2.7 




1 


16 


2.9 


21 


3.8 


80 


3.2 


0.7 




8 


245 


74 


725 


195 


60 


4.9 


1.1 




20 


806 


269 


3442 


967 


40 


7.5 


1.8 




40 


1984 


703 


1.1 X 10" 


2302 


20 


13 


4.7 




60 


3361 


1225 


2.2 X 10" 


6428 


8 


22 


8.0 




80 


4885 


1812 


3.6 X 10" 


1.1 X lO" 


3 


32 


15 




100 


6528 


2451 


5.3 X 10" 


1.5 X 10" 


1 


44 


29 




120 


8274 


3136 


7.2 X 10" 


2.1 X 10" 


0.7 


47 


33 



Table 2. The number fraction ryjv = 100 dm/ J^" ?('^) dm, and the mass 

fraction riM = 100 J"^^ m^{m) dm/Md, M^i = J"^"m^{m)dm, in per cent of 
BDs or main-sequence stars in the mass interval mi to m2, and the stellar con- 
tribution, p^^, to the Oort limit and to the Galactic-disk surface mass-density, 
X"'* = 2 ftp"*, near to the Sun, taking mi = 0. 01 Mq, mu = 120 Mq and the 
Galactic-disk scale-height h = 250 pc (m < 1 Mq iKroupa." Tout fc Gilmore Ill993f ) 
and ft = 90pc (m > 1 Mm. IScalo1ll986l ). Results are shown for the canonical IMF 
(eq. I124p , for the high-mass-star IMF approximately corrected for unresolved com- 
panions (a s = 2 . 7, m > IMq), and for t he present-day mass function (PDMF, 
as = 4.5, IScalo 1 1 19861 : iKroupa. Tout fc"~G ilmore 1993) which describes the dis- 
tribution of stellar masses now populating the Galactic d isk. For gas in th e disk 
^gas ^ ^3 ± sMq/pc^ and remnants E'""" « 3Mq/pc2 (IWeidemannlll990l ). The 
average stellar mass is m = /J^" m£,{m) dm/ £,{m) dm. Nd is the number of 
stars that have to form in a star cluster such that the most massive star in the pop- 
ulation has the mass mmax. The mass of this population is M^i, and the condition 
is (.{fn) dm — 1 with /q™""' dm = A^'d — 1. AMd/Md is the fraction of 

mass lost from the cluster due to stellar evolution, assuming that for m > 8 Mq all 
neutron stars and black holes are kicked o ut due to an asymmetrica l supernova explo- 
sion, but that white dwarfs are retained (|Weidemann et al]|l992l ) and have masses 
mwD = 0.084 mini + 0.444 [Mq]. This is a linear fit to the data in (IWeidemannI 
I2OO0I . their table 3) for progenitor masses 1 < m/MQ < 7 and mwD = 0.5 Mq for 
0.7 < m/MQ < 1. The evolution time for a star of mass mto to reach the turn-off 
age is available in fig. 20 in iKroupa (,2()07ah . 
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expelled by the winds and photo-ionising radiation from massive stars, if they 
are present, being the case in clusters with iV<;few x 100 stars (table [IJ. 
Massive stars (m > 8Mq) are either formed at the cluster ce ntre or get 



there through dynamical mass segregation, i.e. energy equipartition (jBonnell et al 



20071) . The latter process is very rapid (eq.[6l p. [5]) and can occur within 1 Myr 



A cluster core of massive stars is therefore either primordial or forms rapidly 
because of energy equipartition in the cluster, and it is dynamically highly 
unstable decaying within a few icr. core- The ONC, for example, should not 
be hosting a Trapezium as it is extremely volatile. The implication for the 
IMF is that the ONC and other similar clusters and the OB associations 
which stem from them must be very depleted in their massive-star content 



(jPfiamm-Altenburg fc Kroupa II2006I ). 

Important for measuring the IMF are corrections for the typically high 
multiplicity fraction of the very young population. However, these are very 
uncertain because the binary population is in a state of change ffig. [HI below). 

The determination of an IMF relies on the assumption that all stars in a 
very young cluster formed together. However, trapping and focussing of older 
field or OB association stars by the forming cluster has been found to be 
possible (section [TTT]) . 

Thus, be it at the low-mass end or the high-mass end, the stellar mass 
function estimated from very young clusters cannot be the true IMF. Statis- 
tical corrections for the above effects need to be applied and comprehensive 
iV-body modelling is required. 

Old open clusters in which most stars are on or near the main sequence are 
no better stellar samples: They are dynamically highly evolved, since they have 
left their previous concentrated and gas-rich state and so they contain only 



a small fraction of the sta r s originally born in the clus ter (jKroupa fc Boilv 



I2OO2I: IWeidner et al. |[2007i iBaumgardt fc Kroupall2007l) . The binary fraction 



is typically high and comparable to the Galactic field, but does depend on 
the initial density and the age of the cluster, the mass-ratio distribution of 
companions also. So, simple corrections cannot be applied equally to all old 
clusters. The massive stars have died, and secular evolution begins to affect the 
rem aining stellar populatioii (after gas expulsion) through energy equiparti- 
tion. Baumg ardt fc Maki'^ (|2003h have quantified the changes of the MF for 



clusters of various masses and on different Galactic orbits. Near the half-mass 
radius the local MF resembles the global MF in the cluster, but the global MF 
becomes significantly depleted of its low-mass stars already by about 20 per 
cent of the cluster disruption time. 

Given that we are never likely to learn the exact dynamical history of a 
particular cluster, it follows that we can never ascertain the IMF for any indi- 
vidual cluster. This can be summarised concisely with the following theorem: 



Cluster IMF Theorem: The IMF cannot be extracted for any individual 
star cluster. 
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Proof: For clusters younger than about 0.5 Myr star formation has not ceased 

and the IMF is therefore not assembled yet and the clust er cores consisting: of 

mass ive stars have already dynamically ejected members ([Pflamm-Altenburg &: Kroupa 
For clusters with an age between 0.5 and a few Myr the expulsio n 



of residual gas has lead to loss of stars (jKroupa. Aarseth fc Hurley I l200lh . 
Older clusters arc cither still loosing stars due to residual gas expulsion 
or arc evolving secularly through evaporation driven by energy equiparti- 
tion (jBaumgardt fc Makino I l2003h . Furt hermore, the birth sample is likel; 



to be contaminated by captured s tars (jFellhauer. Kroupa fc Evans 1 1200 



I 



Pflamm-Altenburg fc Kroupa Il2007h . There exists thus no time when all stars 
are assembled in an observationally accessible volume (i.e. a star cluster). 

Note that the Cluster IMF Theorem implies that individual clusters 
cannot be used to make deductions on the similarity or not of their IMFs, un- 
less a complete dynamical history of each cluster is available. Notwithstanding 
this pessimistic theorem, it is nevertheless necessary to observe and study star 
clusters of any age. Combined with thorough and realistic TV-body modelling 
the data do lead to essential statistical constraints on the IMF Universality 
Hypothesis (p. [49|) . Such an approach is discussed in the next section. 



The alpha plot 



Scalo I (|l998f ) conveniently summarised a large part of the available observa- 
tional constraint s on the I MF o f resolve d stellar populations with the alpha 
plot, as used by iKroupa I (|200ll . l2002bl ) for explicit tests of the IMF Uni- 
versality Hypothesis (p. SH]) given the Cluster IMF Theorem. One 
example is presented in fig. ^ which demonstrates that the observed scatter 
in a (m) can be readily unde rstood as being due to Poisson uncertainties (see 
also Elmegreen I ll997i Il999l) and dynamical effects, as well as arising from 



biases through unresolved multiple stars. Furthermore, there is no evident 
systematic change of a at a given m with metallicity or density of the star- 
forming cloud. More exotic populations such as the Galactic bulge have also 
been found to have a low -mass MF indistinguishable from the canonical form 
(e.g. lZoccah et al.|[20Q0t ). Thus the IMF Universality Hypothesis cannot 
be falsified for known resolved stellar populations. 



Very ancient and/or metal-poor resolved populations 

Witnesses of the early formation phase of the MW are its globular clusters. 
Such 10^~^Mq clusters formed with individual star- formation rates of 0.1 — 
1 Mq/jv and densities « 5 x 10"^"^ ^Iq/ PC"^. These are relatively high values, 
when compared with the current star-formation activity in the MW disk. For 
example, a 5 x 10^ Mq Galactic cluster forming in 1 Myr corresponds to a star 
formation rate of 0.005 Mq/yt. The alpha plot, however, does not support any 
significant systematic difference between the IMF of stars formed in globular 
clusters and present-day low-mass star formation. For massive stars, it can 
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Fig. 9. The alpha plot. The power-law index, a, is measured over stellar mass-ranges 
and plotted at the mid-point of the respective mass range. The power-law indices 
are measured on the mass function of system masses, where stars not in binaries are 
counted individually. Open circles are the observational constraints for open clusters 
and as s ociati ons for the MW, Large and Small Magellanic clouds collated mostly by 
IScalo I () 19981 ). The open stars (crosses) are theoretical star clusters "observed" in 
the computer at an age of 3 (0) Myr and within a radius of 3.2 pc from the cluster 
centre. The 5 clusters have 3000 stars in 1500 binaries initially and the assumed 
IMF is the canonical one. The theoretical data nicely show a similar spread as the 
observational ones; note the binary-star-induced depression of ai in the mass range 
0.1 — O.5M0. The IMF Universality H ypothesis can t herefore not be discarded 
given the observational data. Models from iKroupa I Hool). 




be argued that the mass in stars more massive than 8 Mq cannot have been 
larger than about half the cluster mass, because otherwise the globular clusters 
would not be as compact as they are today. This constrains the IMF to have 
been close to the canonical IMF (Kroupa 2001). 

A particularly exotic star-formation mode is thought to have occurred in 
dwarf-spheroidal (dSph) satellite galax ies. The MW has abou t 19 such satel- 
lites at distances from 50 to 250 kpc (i Metz fc KroupalboOT ). These objects 
have stellar masses and ages comparable to those of globular clusters but are 
10 — 100 times larger and are thought to have 10 — 1000 times more mass in 
dark matter than in stars. They also show evidence for complex star-formation 
activity and metal-enrichment histories and must have therefore formed un- 
der rather exotic conditions. Nevertheless, the MFs in two of these satellites 
are found to be indistinguishable from those of globular clusters in the mass 
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range 0.5 — 0.9 M(^, thus again showing c onsistency with the canonical IMF 



( Grillmair et al 



The work of 



■■19981: iFeltzing et al.lll999l) . 
Yasui et alj (|2006ll2008[ nias been pushing studies of the IMF 



in young star clusters to the outer, metal-poor regions of the Galactic disk. 
They find the IMF to be indistinguishable, within the uncertainties, with the 
canonical IMF. 



The Galactic bulge and centre 

For low-mass stars the Galactic bulge has been shown to have a MF indis- 
tinguishable from the canonical form ( Zoccali et al. 20001 ). However, abun- 
dance patterns of bulge stars sugg est the IMF to have been top heavy 
(jBallero. Kroupa &: Matteucci Il2007f ). which may be a result of extreme star- 
burst conditions valid in the formation of the bulge ( Zoccali et al. 20061 ). 

Even closer to the Galactic centre, Hertzsprung-Russell-diagram modelling 
of the stellar populatio n within 1 pc of Sgr A * suggests the IMF to have always 



been top-heavy there (jManess et al.l 120071 ). Perhaps this is the long-sought 
after evidence for a variation of the IMF under very extreme conditions, in 
this case a strong tidal field and higher temperatures (but note fig. fTOl below). 



Extreme star bursts 

As noted on p. [2Ql objects with a mass 10^ Mq have an increased M/L 
ratio. If such objects form in 1 to 5Myr, then their star-formation rates, 
SFR^lOM0/yr, and they contain ^10^ O stars packed within a region 
spanning at most a few pc, given their observed present-day mass-radius re- 
lation. Such a star-formation environment is presently outside the reach of 
theoretical investigation. However, it is conceivable that the higher M/L ra- 
tios of such objects may be due to a non-canonical IMF. One possibility is 
that the IMF is bottom heavy as a result of intense photo-destruction of ac - 
cretion envelopes of intermediate to low- mass stars ( Mieske fc Kroupall2008f ). 
Another possibility is that the IMF becomes top-heavy leaving many stellar 
remnants that inflate the M/L ratio ( Dabringhausen fc Kroupa Il2008h . Work 
is in progress to achieve observational constraints on these two possibilities. 



Population III: the primordial IMF 

Most theoretical workers agree that the primordial IMF ought to be top heavy 
because the ambient temperatures were much higher and the lack of metals 
did not a llow gas clouds to cool and to fragment into sufficiently small cores 
( Larson 1 ,1998). The existence of extremely metal-poor low- mass stars with 
chemical peculiarities is interpreted to mean that low-mass stars could form 
under extremely metal-poor conditions, but that their formation was sup- 
pressed in comparison to later star- formation ( Tumlinsonii2007. ) . Modelling of 
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the formation of stellar populations during cosmological structure formation 
suggests that low-mass population III stars should be found within the Galac- 
tic halo if they formed. Th eir absence to-date would imply a primordial IMF 
depleted in low- mass stars ( Brook et al. 2007t ). 

Thus, the last three sub-sections hint at physical environments in which 
the IMF Universality Hypothesis may be violated. 



3.3 Very low-mass stars (VLMSs) and brown dwarfs (BDs) 

The origin of BDs and some VLMSs is being debated fiercely. One camp be- 
lieves these objects to form as stars do, because the star- formation process 
docs not know where the hydrogen burning mass limit is (e.g. Eislof fel fc Steinacker 
I2OQ7). The other camp believes that BDs cannot form exactly like stars 
through continued accretion because the conditi ons required for this to 



occur in molecular clouds a re far too rare (e.g. iReipurth &: Clarkd 12001 



Goodwin fc Whitworthll2007 ) 



If BDs and VLMSs form like stars then they should follow the same pair- 
ing rules. In particular, BDs and G dwarfs would pair in the same manner, 
i.e. according to th e same mathematical rules, as M dwarfs and G dwarfs. 



Kroupa et al.l (|20031 ) tested this hypothesis by constructing iV-body models of 



Taurus- Auriga-like groups and Orion-Nebula-like clusters finding that it leads 
to far too many star-BD and BD-BD binaries with the wrong semi-major axis 
distri bution. Instead, star-BD binaries are very rare (jGrether fc Lineweaver 
while BD-BD binaries are rarer than stellar binaries (BDs have a 15 per 



cent binary fraction as opposed to 50 per cent for stars), and BDs have a semi- 
major axis distribution significantly narrower than that of star-star binaries. 
The hypothesis of a star-like origin of BDs must therefore be discarded. BDs 
and some VLMSs form a separate population, which is however linked to that 
of the stars. 

iThies fc Krounal (120071 ) re-address this problem with a detailed analysis of 
the imderlying MF of stars and BDs given observed MFs of four populations, 
Taurus, Trapezium, IC348 and the Pleiades. By correcting for unresolved bi- 
naries in all four populations, by taking into account the different pairing rules 
of stellar and VLMS and BD binaries, a significant discontinuity of the MF 
emerges. BDs and VLMSs therefore form a truly separate population from 
that of the stars and can be described by a single power-law MF (eq. I129|) 
which implies that about one BD forms per 5 stars in all four populations. 

This strong correlation between the number of stars and BDs, and the 
similarity of the BD MF in the four populations implies that the formation of 
BDs is closely related to the formation of stars. Indeed, the truncation of the 
binary binding energy distribution of BDs at a high energy suggests that en- 
ergetic processes must be operating in the production of BDs, as discussed in 
iThies fc Kroupa ( 2007). T wo such possible mechanis ms are embryo ejection 
( Reipurth fc Clarkd 2001 ) and disk fragmentation ( Goodwin fc Whit worth 
20071 ). 
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3.4 Composite populations: the IGIMF 

The vast majority of all stars form in embedded clusters and so the correct 
way to proceed to calculate a galaxy-wide stellar IMF is to add-up all the 
IMFs of all star-clust ers born in one st ar- formation epoch. Such epochs may be 
identified with the Zoccali et al. ( 20061 ) star-burst events creating the Galactic 



bulge. In disk galaxies they may be related to the time-scale of transforming 
the inter-stellar matter to star clusters along spiral arms. Addition of the 
clusters bor n in one epoch gives the in tegrated galactic initial mass function, 
the IGIMF (iKroupa fc Weidner Il2003h . 



IGIMF Definition: The IGIMF is the IMF of a composite population which 
is the integral over a complete ensemble of simple stellar populations. 



Note that a simple population has a mono-metallicity and a mono-age distri- 
bution and is therefore a star cluster. Age and metallicity distributions emerge 
for massive populations with M^i ^ 10^ Mq (e.g. oj Gen) indicating that such 
objects, which also have relaxation times comparabe to or longer than a Hub- 
ble time, are not simple (section II. 4p . A complete ensemble is a statistically 
complete representation of the initial cluster mass function (ICMF) in the 
sense that the actual mass function of A'ci clusters lies within the expected 
statistical variation of the IGMF. 



IGIMF Theorem: The IGIMF is steeper than the canonical IMF if the IMF 
Universality Hypothesis holds. 



Proof: IWeidner k Kroupa I (|2006l ) calculate that the IGIMF is steeper than 



the canonical IMF for m > 1 M© if the IMF Universality Hypothesis 
holds. The steepening becomes negligible if the power-law mass function of 
embedded star clusters, 

^oci(Mcci) oc mJ (127) 

is flatter than (3 = 1.8. 

It may be argued that IGIMF=IMF (e.g. Elmegreen 2006) because when 
a star cluster is born, its stars are randomly sampled from the IMF up to 
the most massiv e star possible . On the oth er hand, the physically- motivated 
ansatz by Weidne r fc Kroupa ( 20051 [2006h of taking the mass of a cluster as 



the constraint and of including the observed correlation between the maxi- 
mal star mass and the cluster mass (fig. [1]), yields an IGIMF which is equal 
to the canonical IMF for m^l.SM© but which is systematically steeper 
above this mass. By incorporating the observed maximal-cluster-mass vs star- 
formation rate of gala xies, Mpc]. rria.^ = M c;r.].ma.x(SFR), for the youngest clus- 



ters (Weidner, Kroup a fc LarserLll2004( l. it follows for Mq that low- 



surface-brightness (LSB) galaxies ought to have very steep IGIMFs, while 
normal or L* galaxies have Scalo-type IGIMFs, i.e. aiciMF = aiviwdisk > ^2 
(section l3.ip follows naturally. This systematic shift of aiciMF (m^l.SM©) 
with galaxy type implies that less-massive galaxies have a significantly sup- 
pressed supernova II rate per low-mass star. They also show a slower chem- 
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ical enrichment such tha t the observed metallicity-galaxy-m ass relation can 
be nicely accounted for ( Koeppen. Weidner &: Kroupa 112007 ). Another very 
important implication is that the SFR-Ha-luminosity relation for galaxies 
flattens such that the SFR becomes higher by up to three orders of magni- 
tude for dwarf gala xies than the value calculated from the stand ard (linear) 
Kennicutt relation ('Pflam m-Altenburg. Weidner &: Kroupal 2007 ). 



St rikingly, the IGIMF variation has now been direct ly measured bv | Hoversten &: Glazebrook 
( 2008h using galaxies in the Sloan Digital Sky Survev. [Lee et al" have 
inde ed found LSBs to hav e bottom- heavy IMFs, while IPortinari et al.l ( 2004 ) 
and Romano et al. I liooi) find the MW disk to have a steeper- t han S alpeter 



IMF for massive stars which is, in comparison with iLee et al. I (|2004| ). much 



flatter than the IMF of LSBs, as required by the IGIMF Theorem. 



3.5 Origin of the IMF: theory vs observations 



General physical concepts such as coalescence of proto-stellar cores, mass- 
dependent focussing of gas accretion onto proto-stars, stellar feedback, and 
fragmentation of molecular clouds lead to predictions of systematic varia- 
tions of the IMF with changes of the physical conditions of sta r forma- 



tion (iMurrav fc LinI llQQd lElmegreenI |2004 iTillev fc Pudritd l2005i but see 

i 



Casuso fc Beckmanll2007l for a simple cloud coagulation/dispersal model lead- 
ing to an invariant mass distribution). Thus, the thermal Jeans mass of a 
molecular cloud decreases with temperature and increasing density, implying 
that for higher metallicity {— stronger cooling) and density the I MF should 
shift on average to smaller stellar masses (e.g. Lar son 1998; Bonnell et al 



20071 ). The entirely different notion that stars regulate their own masses 



through a balance between feedback and accretion also implies smaller stellar 
masses for higher metallicity due in part to more dust and thus more efficient 
radiation pressure on the gas through the dust grains. Also, a higher metal- 
licity allows more efficient cooling and thus a lower gas temperature, a lower 
sound speed and therefor e a lower accretion rate (Adams & Fatuzzo 19961; 



Ed ams fc Laughlin I Il996l ). As discussed above, a systematic IMF variation 
with physical conditions has not been detected. Thus, theoretical reasoning, 
even at its most elementary level, fails to account fo r the observat i ons. 

(l2007f ) re- 



Klessen et al 



A dramatic case in point has emerged recently: 
port state-of-the art calculations of star-formation under physical conditions 
as found in molecular clouds near the Sun and they are able to reproduce 
the canonical IMF. Applying the same computational technology to the con- 
ditions near the Galactic centre they obtain a theoretical IMF in agree- 
ment with the previously reported apparent d e cline of the stellar MF in the 
Arches cluster below about 6Mq. iKim et aL ( 2006h pubhshed their obser- 
vatio ns of the Arches clu ster on the astro-physics preprint archive shortly 
after iKlessen et al. I (l2007l) and performed the necessary state-of-the art A^- 
body calculations of the dynamical evolution of this young cluster, revising 
our knowledge significantly. In contradiction to the theoretical prediction they 
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find that the MF continues to increase down to their 50 per cent complete- 
ness limit (1.3 Mq) with a power-law exponent only slightly shallower than 
the canonical Massey/Salpeter value once mass-segregation is corrected for. 
This situation is demonstrated in fig. [TOl 



□ 




log (M/Mg) 

Fi g. 10. The observed mass function of the Arches cluster near the Galactic centre 
bv lKim et al.l (|2006i ) shown as the thin histogramme is confro nted with the theoret - 
ical MF for this object calculated with the SPH technique bv lKlessen et al.l (|2007h . 
marked as the hatched histogramme. The latter has a down-turn (bold steps near 
10'''^) incompatible to the observations therewith ruling out a theoretical under- 
standing of the stellar mass spectrum (one counter-example suffices to bring-down 
a theory). One possible reason for the theor etical failure may be the assumed tur- 
bulence driving. For details on the figure see iKim et al.l l|2006h . 



It therefore emerges that there does not seem to exist any solid theoretical 
understanding of the IMF. 

Observations of cloud cores appear to suggest t hat the cano nical IMF is 
frozen-in already at the pre - stella r clo ud-core level ( Motte et al][T99 8. 200j). 
Nutter fc Ward-ThompsonI (|2007l ) and lAlves et al.l (|2007t) find, however, the 



pre-stellar cloud cores to be distributed according to the same shape as the 
canonical IMF, but shifted to larger masses by a factor of about three or 
more. This is taken to perhaps mean a star-formation efficiency per star of 
30 per cent or less independently of stellar mass. The interpretation of such 
observation s in view of mu l tiple star for mation in e ach cl oud-core is being 
studied bv I Goodwin et aTl ( 2008f ). while Krumhol3 ( 2007 ) outlines current 
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theoretical understanding of how massive stars form out of massive pre-stellar 
cores. 



3.6 Conclusions: IMF 

The IMF Universality Hypothesis, the Cluster IMF Theorem and 
the IGIMF Theorem have been stated. Furthermore, 

1. The stellar luminosity function has a pronounced maximum at My « 
12, Ml K, 9 which is universal and well understood as a result of stellar 
physics. Thus by counting stars on the sky we can look into their interiors. 

2. Unresolved multiple systems must be accounted for when the MFs of 
different stellar populations are compared. 

3. BDs and some VLMSs form a separate population which correlates with 
the stellar content; there is a discontinuity in the MF near the star/BD 
mass transition. 

4. The canonical IMF (eq. I124p fits the solar-neighbourhood star counts and 
all resolved stellar populations available to-date. Recent data at the Galac- 
tic centre suggest a top-heavy IMF, perhaps hinting at a possible variation 
with conditions (tidal shear, temperature). 

5. Simple stellar populations are found in individual star clusters with Afd 
< 10'^ Mq. These have the canonical IMF. 

6. Composite populations describe entire galaxies. They are a result of many 
epochs of star-cluster formation and are described by the IGIMF Theo- 
rem. 

7. The IGIMF above ~ 1 Mq is steep for LSB galaxies, flattening to the 
Scalo slope (aiciMF ~ 2.7) for L* disk galaxies. This is nicely consistent 
with the IMF Universality Hypothesis in the context of the IGIMF 

THEOREM. 

8. Therefore, the IMF Universality Hypothesis can not be excluded 
despite the cluster IMF Theorem for conditions p^lO^ stars/pc^, 
Z<; 0.002 and non-extreme tidal fields. 

9. Modern star-formation computations and elementary theory give wrong 
results concerning the variation and shape of the stellar IMF, as well as 



the stellar multiplicity (jGoodwin fc Kroupa Il2005f ) 



10. The stellar IMF appears to be frozen-in at the pre-stellar cloud-core stage 
therewith probably being a result of the processes leading to the formation 
of self-gravitating molecular clouds. 



3.7 Discretisation 

As discussed above, a theoretically-motivated form of the IMF which passes 
observational tests does not exist. Star-formation theory gets the rough shape 
of the IMF right (there are fewer massive stars than low-mass stars), but other 
than this, fails to make any reliable predictions whatsoever as to how the IMF 
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should look like in detail under different physical conditions. In particular, the 
overall change of the IMF with metallicity or density or temperature predicted 
by theory is not evident. An empirical multi-power-law form description of the 
IMF is therefore perfectly adequate, and has important advantages over other 
formulations. A general formulation of the stellar IMF in terms of multiple 
power-law segments is as follows: 



f (m) = fc< 



-ao 
-ai 



\rr 

( ™o ^ ^ ( J2l\ 
\ mu J \mo J 

fe) te) 



miow < m < mn 
mn < m < mo 
mo < m < mi 



(128) 



where mmax < Timax* ~ 150 M0 depends on the stellar mass of the embedded 
cluster (fig. [Ij. The empirically determined stellar IMF is a two-part- form 
(eq. I124p . with a third power-law for BDs, whereby BDs and VLMSs form a 
separate population from that of the stars fp. [55]) . 



CBDOcm~"«, ao~0.3, (129) 
(jMartfn et al.ll200oHchabrier |[200l iMoraux et al.ll2004l )) and 



Cbd (0.075 Mq) w 0.25 C(0.075Me 



(|Thies fc Kroupall2007h where ^ is the canonical stellar IMF feq. \T^ . This 
implies that about one BD forms per 5 stars. 

One advantage of the power-law formulation becomes immediately obvious 
by realising that analytical generating functions, and other quantities, can be 
derived readily. Another important advantage is that by using a multi-power- 
law form, different parts of the IMF can be varied in numerical experiments 
without affecting the other parts. A practical num e rical formulation of the 
IMF is prescribed in Pflamm-Altenburg fc Kroupa I ( 20061 ). Thus, for exam- 
ple, the canonical two-part power-law IMF can be changed by adding a third 
power-law above 1 Mq and making the IMF top-heavy {am>i Mq < 0^2) with- 
out affecting the shape of the late-type stellar luminosi ty function as evident 
in fig. E The KTG93 (|Kroupa. Tout fc Gilmore1ll993f ) IMF is such a three- 
part power-law form relevant for describing the overall young population in 



the M ilky Way disk, which is top-light {am>iAiQ > a-? . iKroupa fc Weidner 
20031) . 

A log-normal formulation does not offer these advantage, and requires 
power-law tails above about IM©, and for brown dwarfs, for consistency 
with the observational constraints discussed above. However, while not as 
mathematically c onvenient, the p opular Chabrier log-normal plus power-law 
IMF (table 1 in IChabrier" 120031 ) formulation leads to an indistinguishable 



stellar mass distribution to the two-part power-law IMF (fig. [TT|) . Various 




Fig. 11. Comparison between the popular Chabrier IMF (loK-norma l plu s 
power-law extension above 1 Mq : dashed curve, table 1 in Ichabrier I |2003| ) 
with the canonical two-part power- law IMF (solid line, eq. I124p . From 
IPabringhausen Hilker fc Kroupa I (|2008h . 



analytical forms for the IMF are compiled in table 3 of iKroupal (|2007a[ ). 

A generating function for the two-part power-law form of the canonical 
IMF, eg 11241 can be written down by following the steps taken in section [ 
The corresponding probability density is 



Pi - m-^i , 0.08 < m < 0.5 Mq 

P2 = fcp,2 m "^ 0.5 < TO < TOmax, 



where kp^i are normalisation constants ensuring continuity at 0.5 Mq and 

pi dm + P2 dm — 1, (131) 

0.08 Jo. 5 

whereby TO,„ax follows from fig. [T] Defining 

<.0.5 

X[= pi{m)dm, (132) 
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it follows that 



or 



Xi{rn) = / pi{m) dm, ifm<0.5M, 
Jo. 08 

r 

X^im) =X[ + 

Jo. 



P2{m) dm, ifTO>O.5M0. 



(133) 



(134) 



The generating function for stellar masses follows by inverting the above two 
equations Xi{m). The procedure is then to choose a random variate X G [0, 1] 
and to select the generating function m{Xi ~ X) if < X < Xi, or m{X2 = 
X) if Xi < X < 1. 

This algorithm is readily generalised to any number of power-law segments 
(eq. I128p , such as including a third segm ent for brown dwarfs an d allowing the 
IMF to be discontinuous near 0.08 (jThies fc KroupallioOTh . Such a form 



has been incorporated into Aarseth's Nbody4/6/7 programmes, but hitherto 
without the discontinuity. However, Jan Pflamm-Altenburg developed a more 
powerful and general method of generating ste llar masses (or any other quan- 
tities) given an arbitrary distribution function (jPflamm-Altenburg fc Kroupa 
l2OO6.Pl . 



The C-language software package "libimf", can be downloaded 
jhttp : //www, astro .uni-bonn. de/~webaiub/english/downloads .phpl . 



at 
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4 The initial binary population 
4.1 Introduction 

It has already been demonstrated that corrections for unresolved multiple 
stars are of much importance for correctly deriving the shape of the stellar MF 
given an observed LF (fig. [5]). Binary stars are also of significant importance 
for the dynamics of star clusters, because a binary has intrinsic dynamical 
degrees of freedom which a single star does not possess. A binary can therefore 
exchange energy and angular momentum with the cluster. Indeed, binaries 
are very significant energy sources, as for example, a binary composed of two 
1 Mq main sequence stars and with a semi- major axis of 0.1 AU has a binding 
energy comparable to that of a 1000 Af© cluster of size 1 pc. Such a binary 
can interact with cluster-field stars accelerating them to higher velocities and 
thereby heating the cluster. 

The dynamical properties describing a multiple system are: 

• the period P (in days throughout this text) or semi-major axis a, 

• the system mass Wgys = mi + m2, 

• the mass ratio q = ^ < 1, where TOi,m2 are, respectively, the primary 
and secondary-star masses, 

• the eccentricity e = (rapo - r-pori)/(rapo + »^pori), where rapo,rpori are, re- 
spectively, the apo-centric and peri-centric distances. 

Given a snapshot of a binary, the above quantities can be computed from 
the relative position, rid, and velocity, v^ci, vectors and the masses of the two 
companion stars by first calculating the binding energy, 

1 2 Gmim2 Gniim2 

2 Ti-ci 2 a 

where = mi m2 / {mi + TO2) is the reduced mass. From Kepler's third law 
we have 

TUsys = =^ P = Py, X 365.25 days, (136) 
yr 

where Pyr is the period in years and cau is the semi-major axis in AU. Finally, 
the instantaneous eccentricity can be calculated using 



2 ^ (r-rcl ■ ^'rel)^ 



(137) 



which can be derived from the orbital angular momentum, 

i = /Xt>rel X rrel, (138) 



with 
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L 

'^^sys 

The relative equation of motion is 



— a(l-e2) 



17111712- (139) 



= -G^r,.ei + apert(0, (140) 

rcl 



where apert(i) is the time-dependent perturbation from other cluster members. 
It follows that the orbital elements of a binary in a cluster are functions of time, 
P — P{t) and e — e{t). Also, q = q{t) during strong encounters when partners 
are exchanged. Since most stars form in embedded clusters it follows that the 
binary-star properties of a given population cannot be taken to represent the 
initial or primordial values. 

The following theorem can therewith be stated: 



Dynamical Population Synthesis Theorem: If initial binary populations 
are invariant, then a dynamical birth configuration of a stellar population can 
be inferred from its observed binary population. This birth configuration is not 
unique, however, but defines a class of dynamically equivalent solutions. 



The proof is simple: Set-up initially identical binary populations in clusters 
with different radii and masses, and calculate the dynamical evolution with an 
TV— body programme. For a given snapshot of a p opulation, ther e is a scalable 
starting configuration in terms of size and mass (lKroupalll9953 lcih. 



Binaries can absorb energy and thus cool a cluster. They can also heat 
a cluster. There are two extreme regimes that can be understood with a 
Gedanken experiment. Define 

i?k = (1/2) ma^ w (1/iV) X kinetic energy of cluster. U^-Lj 

Soft binaries have i?bin ^ -E-k, while hard binaries have £'bin S> -Bk- A useful 
equation in this context is the relation between the orbital period and circular 
velocity of the reduced particle F^: 

/P[days] = 6.986 + ZmsysiM©] - 3 /uorb[km/s], (142) 

Consider now the case of a soft binary, i.e. a reduced-mass particle with 
Vorh ^ cr. By the principle of energy equipartition, Worb c (eq. [5]) as time 
progresses. This implies a '\,P t- A hard binary has Uorb 3> cr. Invoking 
energy equipartition, it follows that Uorb i, i-e. a |,P |. Furthermore, the 
amount of energy needed to "ionise" a soft binary is negligible compared to 
the amount of energy required to "ionise" a hard binary. And, the cross section 
for suffering an encounter scales with the semi-major axis implying that a soft 
binary becomes ever more likely to suffer an additional encounter as its semi- 
major axis increases. Therefore, it is much more probable for soft binaries 



Throughout this text: Ix = logj^Q(2;). 
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to be disrupted rapidly, than for hard binary to do so. Thus follows ( Heggid 
197 5 iHilldllQTSh 



Heggie-Hills Law: Soft binaries soften and cool a cluster. Hard binaries 
harden and heat a cluster. 



Numerical scattering experiments by 'Hilld ( 1975) have shown that harden- 



ing of binaries often involves partner exchanges. iHeggie derived the 

above law analytically. Binaries in the energy range 10~^ £^k ~ -^bin S5 10^ E'k, 
33~^ cr ^ Worb ~ 33 cr cannot be treated analytically due to the complex reso- 
nances that are created between the binary and incoming star or binary. It is 
these binaries that may be important for the early cluster evolution, depend- 
ing on its velocity dispersion, a = (t(Moci). Cooling of a cluster is energetically 
not significant but has been seen for the first time by Kroupa et al. (_1999). 

Fig. [m visualises the broad evolution of the initial period distribution in a 
star cluster. At any time, binaries near the hard/soft boundary, with energies 
E'bin ~ a-nd periods P ~ Pth(worb = o") (eq. El) denoting the thermal 
period^ are most active in the energy exchange between the cluster field and 
the binary population. The cluster expands as a result of binary heating and 
mass segregation, and the hard/soft boundary, Pth, shifts to longer periods. 
Meanwhile, binaries with P > Pth continue to be disrupted while Pth keeps 
shifting to longer periods. This process ends when 

Pth^Pcut, (143) 

which is the cutoff or maximum period in the surviving period distribution. 
At this critical time, ^i, further cluster expansion is slowed because the pop- 
ulation of heating sources, i.e. the binaries with P « Pth, is significantly 
reduced. The details strongly depend on the initial value of Pth which deter- 
mines the amount of binding energy in soft binaries which can cool the cluster 
if significant enough. 

After the critical time, tt, the expanded cluster reaches a temporary state 
of thermal equilibrium with the remaining binary population. Further evolu- 
tion of the binary population occurs with a significantly reduced rate deter- 
mined by the velocity dispersion in the cluster, the cross section given by the 
semi-major axis of the binaries, their number density and that of single stars 
in the cluster. The evolution of the binary-star population during this slow 
phase will usually involve partner exchanges and unstable but also long-lived 
hierarchical systems. The IMF is critically important for this stage, as the ini- 
tial number of massive stars determines the cluster density at t ^ 5 Myr owing 
to mass loss from evolving stars. Further binary depletion will occur once the 
cluster goes into core-collapse and the kinetic energy in the core rises. 
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At t=tt 




^ [d] 



for a3jj=3.5 km/s 
at f=0 

Fig. 12. Illustration of the evolution of the distribution of binary star periods in a 
cluster {IP — logigP). A binary has orbital period Pth when usd (= f) equals its 
circular orbital velocity (eq. I142p . The initial or birth distribution (eq. I165P evolves 
to the form seen at time t > tt- 



Frequency of binaries and higher-order multiples 



The emphasis here is on late-type binary stars because higher-order multi- 
ples are rare as shown by obs ervation. The iiiforma tion on the multiplicity of 
massive stars is very limited (jGoodwin et al.ll2007l ). 

Defining, respectively, the number of single stars, binaries, triples, quadru- 
ples, etc., by the numbers 



(iVsing : Wbin : Nt, 



ip 



quad 



= {S : B : T : Q : ...) (144) 



the multiplicity fraction can be defined, 

Nn,uit B + T+Q 



fmult — 



sys 



S+B+T+Q 



(145) 



and the binary fraction is 



Dense Stellar Populations: Initial Conditions 
B 



67 



/bin 



TVs, 



(146) 



In the Galactic field iDuauennov fc Mavoi ( 199l[) derive from a decade- 
long surv ey for G-dwarf primary s tars, ^A^muit = ( 57 : 38 : 4 : 1 ) and for 
M-dwarfs lrTscher fc Marcvl (|l992h find ^TVmuit = ( 58 : 33 : 7 : 1 ). Thus, 



^fmnlt = 0.43; 
^/muit = 0.41; 



"^/bin = 0.38 
^/bin = 0.33 



It follows that most "stars" are indeed binaries. 
After correcting for incompleteness, 



0.53 ±0.08, 



7bin = 0.45 ±0.07, 



M , 



(150) 
(151) 

Vbin = 0.42 ± 0.09, (152) 

where the K-dwarf data have been published bv lMavor et al. ( 19921 ). It follows 
that 

^/bin ~^ /bin /bin ~ 0.5 « /tot, (153) 

in the Galactic field, perhaps with a slight decrease to wards lower masses. In 



contrast, for brown dwarfs, ^°/bin ~ 0.15 /bin (|Thies fc Kroupall2007 . 

and references therein). 

An interestin g problem arise s when one considers that for 1 Myr old stars, 
/xTauri ~ 1 (s-g. Ducheud 19991 ). Given the above information, the following 
theorem can be stated: 



Binary-Star Theorem: The vast majority of stars form in binary systems. 



Proof: If a substantial fraction of stars were to form in higher-order multiple 
systems, or as small— systems, then the typical properties of these at birth 
imply their decay within typically 10"* to 10^ yr leaving a predominantly single- 
stellar population. Ho wever, the majority of 10^ M yr old stars are observed to 
be in binary systems ( Goodwin fc Kroupa 2005f ). 

Higher-order multiple systems do exist and can only be hierarchical to 
guarantee stability. Hierarchical systems are multiple stars which are stable 
over many orbital times, and are typically tight binaries orbited by outer 
tertiary companions , or tw o tight binaries in orbit about each other, etc. 
Eggleton fc Kiseleva (1993) discuss stability issues in more detail, but it suf- 
fices here to state that the outer and inner semi-major axes should typically 
have a ratio of about four for stability for comparable system masses. If the 
stability criterion is not fulfilled, then higher-order multiple systems typically 
decay on a few crossing times by ejecting members until a stable or long-lived 
configuration is found. Most often this is a hardened binary. 

Star cluster remnants (or dead star clusters) may be the origin of most 
hierarchical, higher-order multiple stellar systems (p. I20|). 
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4.2 The initial binary population late-type stars 



The initial binary population is described by distribution functions that are as 
fundamental for a stellar population as the IMF. There are four distribution 
functions defining the initial dynamical state of a population: 

1. The IMF, f(m), 

2. the distribution of periods (or semi-major axis), df — fp[lP) dlP 

3. the distribution of mass-ratios, df — fq{q) dq, 

4. the distribution of eccentricities, df = fe{e)de, 

where df is the fraction of systems with a parameter in the vicinity of the 
given value. Thus, for example, fip{lP = 4.5) = 0.11, i.e. of all G-dwarfs on 
the sky, 11 per cent have a companion with a period in the range of 4 to 5 days 

These distribution functions have been measured for late-t ype stars in the 

Galac tic fie ld and in star-form i ng reg ions ffig. [T3|) . According to lDuquennov &: Mayor 
(1991) and iFischer fc Marcvl ( 19921 ) both G-dwarf and M-dwarf binary sys- 
tems in the Galactic field have period distribution functions that are well 
described by log-normal functions. 



fp{lP) = fu 



1 



(TIP V2 vr 



(154) 



with IP « 4.8 and aip « 2.3, and J^,; ,p fip{lP) dlP = /tot ~ 0.5. K-dwarfs 
appear to have an indistinguishable period distribution. 

From fig. [13] it follows that the pre-main sequence binar y fraction is larger 
than that of main-sequence stars (see also Duchenelll999f) . Is this an evolu- 
tionary effect? 

Further, Duquennov fc Mavod (1991) derived the mass-ratio and eccen- 
tricity distributions for G-dwarfs in the Galactic field. The mass-ratio dis- 
tribution of G-dwarf primaries is not consistent with random sampling from 
the canonical IMF (eq. I124p. as the nu mber of observed low-mass compan- 
ions is underrepresented (|Kroupalll995cl ). In contrast, the pre-main sequence 



mass-ratio distribution is consistent, within the uncertainties, with random 



sampling from the canonical IMF for g^O.2 (|Woitas et al.ll2001r ). The eccen- 
tricity distribution of Galactic-field G-dwarfs is found to be thermal for IP ^ 3, 
while it is bell shaped with a maximum near e — 0.25 for IP ^ 3. Not much is 
known about the eccentricity distribution of pre-main sequence binaries, but 
numerical experiments show that fe does n ot evolve much in dense clusters, 
i.e. the thermal distribution must be initial ( Kroupa 1995dl ). 
The thermal eccentricity distribution, 



/e(e) = 2e. 



(155) 



follows from a uniform binding-energy distribution (i.e., all energies are 
equally populated), as follows. The orbital angular momentum of a binary 
is 
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Fig. 13. Measured period-distributio n functions for Galactic- field G-dwarfs (his - 
togramme, iDuguennov fc Mayor 1991), K-dwar fs (open circles, iMavor et al.lll9Sj3 ) 
and M-dwarfs (asterisks. iFischer fc Marcvl[l993 ). The about IMyr old T Tauri bi- 
nary data (open squares, partially from the T aurus- Auriga stellar Rroups) are a 
compilation from various sources (see fig. 10 in lKroupa. Aarseth fc Hurley I I2OOII ). 
In all cases, the areas under the distributions is /tot- 



G G mi 1712 



^sys 



2i?bii 



from which follows 



Differentiation leads to 
de 



1 — 2 i?bin L 



G^ (mi TO2) 



(1 — e^) (mi 

mgys 
G^ (mi m2)^ , 



(156) 



(157) 



(158) 
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The number of binaries with eccentricities in the range e,e + de is the same 
number of binaries with binding-energy in the range -Bbinj £'bin + c'-^'bin (the 
same sample of binaries), 

/(e) de = /(£;bi„) dEun = /(^bin) []-^ede, (159) 

where the square brackets are from the previous equation. But 

1 

f{e)de = l, (160) 



I.e. 



r.1 

1 - /(^bin) / ede. (161) 



So 



/(-Ebin) []"^ = 2 = const =^ f{e)de^2ede. (162) 

Thus, /(e) = 2 e is a thermalised distribution: all energies are equally occupied 
{f{Eun) = const). 

iV— body experiments have demonstrated that the period distribution 
function must span the observed range of periods at birth, as dynamical 
encounters in dense clus ters cannot widen an initially narrow distribution 
(|Kroupa k BurkertllioOll ). 

There are thus three discrepancies between main-sequence and pre-main 
sequence late-type stellar binaries: 

1. the binary fraction is higher for the latter, 

2. the period distribution function is different, and 

3. the mass-ratio distribution is consistent with random paring for the latter, 
while it is deficient in low-mass companions in the former, for G-dwarf 
primaries. 

Can these be unified, i.e. are there unique initial //p, /q, /e consistent with the 
pre-main sequence data that can be evolved to the observed main-sequence 
distributions? 

This question can be solved by framing the following ansatz: Assume the 
orbital-parameter distribution function for binaries with primaries of mass mi 
can be separated, 

V{lP,e,q : mi) = fipf^U (163) 

The stellar- dynamical operator, Sl^'^° '^, can now be introduced such that the 
initial distribution function is transformed to the final (Galactic-field) one, 

DnnilP,e,q : mi) = = [V,,,{lP,e,q : mi)] . (164) 

This operator provides a dynamical environment equivalent to that of a star 
cluster with N stars and a half-m ass radius rn..»s (s ee al so the Dynamica l Popu- 
lation Synthesis Theorem, p.[64|). Kroupa (1995c) and Kroupa ( 1995dl ) indeed 
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prove this to be the case for a cluster N = 200 binaries and rg.s — 0.77 pc and 
derives the initial distribution function, for late- type binary systems such 
that it fulfills the above requirement and also has a simple generating function 
(see below). Noteworthy is that such a cluster is very similar to the typical 
cluster from which most field stars probably originate. The full solution for 
fi, such that the Galactic-field is re-produced from forming and dissolving 
star clusters, requires full-scale inverse dynamical population synthesis for the 
Galactic field. 

Thus, by the Dynamical Population Synthesis Theorem (p. [Ml), the 
above ansatz with fi^'^"-^ leads to one solution to inverse dynamical popula- 
tion synthesis (the 200 binary, ro.5 — 0.8 pc cluster, fig. [Mj i.e. most stars in 
the Galactic field stem from clusters dynamically similar to this one) , provided 
the birth (or primordial) distribution functions for /P, e, q are as follows: 

IP —IP ■ 



This distribution function has a generating function (section [ 



The solution obtained by iKroupa (1995d) has 



^ = 2.5, 5 = 45, /P,ni„ = l, (167) 

such that ^-Pmax = 8.43 since //^f'"'"' jip dlP — /tot = 1 is a requirement 
for stars at birth. Intriguingly, similar distributions can be arrived at semi- 
empirically assu ming an isol ated formation of binary stars in a turbulent 
molecular cloud ( Fisher 20041 ). 



The birth-eccentricity distribution is thermal (eq. I155P while the birth 
mass-ratio distribution is generated from random pairing from the canonical 
IMF. However, in order to re-produce (1) the observed data in the eccentricity- 
period diagramme, (2) the observed eccentricity distribution and (3) the ob- 
served mass-ratio distribution for short-period {IP ^3) systems, a correlation 
of the parameters needs to be introduced through eigenevolution. Eigenevolu- 
tion is the sum of all dissipative physical processes that transfer mass, energy 
and angular momentum between the companions when they are still very 
young and accreting. 

A formulation which is quite successful in re-producing the overall observed 
correlations between /P, e, q for short-peri od systems has been derived from 
tidal-circularisation theory ( Kroupal 1995dl ). Most-effective orbital dissipation 



occurs when the binary is at its peri-astron, 

r-pcri = (l-e)Pyl (mi-fma)^ (168) 

where Pyr — P/365.25 is the period in years. Let the binary be born 
with eccentricity ehirt.hi th en the system evolves approximately according to 



(Goldman fc Mazeh 



1994 ) 
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Fig. 14. Evolution of the total binary fraction for stellar mass 0.1 < mi/MQ < 
1.1, i = 1, 2 with time for the four star-cluster models initially with A*' = 200 binaries 
computed by iKroupal (|l995cl l in the search of the existence of an ^2""°^'^ (-R0.5 is 
the initial half-mass radius of the clusters, denoted in this text as ro.s). Note that 
the ro.5 = 0.8 pc cluster yields the correct /tot ~ 0.5 for the Galactic field. The 
period-distribution function and the mass-ratio distribution function that emerge 
from this cluster also fit to the observed Galactic-field distribution. Some binary 
stars form by three-body encounters in clusters that initially consist only of single 
stars, and the proportion of such binaries is shown for the single star clusters (with 
initially A'^ = 400 stars). Such dynamically formed binaries are very rare and so /tot 
remains negligible. 
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1 de 

- ^ = -p' =^ logioeiii = -p + logioebirth, (169) 

where 1/p' is the tidal-circularisation time-scale, ein is the initial eccentricity, 
and 



p'dt=[^\ , (170) 

where Rq is the Solar radius in AU, A, x are tidal-circularisation parame- 
ters and Tpori (in AU) is assumed to be constant because the dissipational 
force only acts tangentially at peri-astron. Note that a large A implies that 
tidal-dissipation is effective for large separations of the companions (e.g. they 
are puffed-up pre-main sequence structures), and a small x implies the dis- 
sipation is soft, i.e. weakly varying with the separation of the companions. 
In this integral. At ^ 10^ yr is the time-scale within which pre-main-sequence 
eigenevolution completes. The initial period becomes, from eq. 11681 



Pin ^ PbUth ('!^^^^) " ( 1^11^] ' , (171) 



Kroupal (|l995dh assumed the companions merge if a-m < 10 Rq : mi 



-TO2 



m. 

In order to re-produce the observed mass-ratio distribution given random 
pairing at birth, and to also re-produ ce the fa c t that short-period binaries 
tend to have similar-mass companions, Kroupal ( 1995d) implemented a feed- 



ing algorithm, according to which the secondary star accretes high-angular- 
momentum gas from the circum-binary accretion material, such that its mass 
increases while the primary mass remains constant. Thus, after generating 
the two birth-companion masses randomly from the canonical IMF, the ini- 
tial mass-ratio becomes 



gin = %irth + (1 - gbirth) P* , (172) 

where 

The above is a very simple algorithm which nevertheless re-produces the 
essence of orbital dissipation such that the correlations between the orbital pa- 
rameters for short-period systems are well accounted for. The best parameters 
for the evolution "birth — > initial" are 



A = 28, x = 0.75. (174) 

Fig- ESI shows an example of the overall model in terms of the eccentricity- 
period diagramme. Fig. [TBI and fig [T71 demonstrate that it nicely accounts for 
the period and mass-ratio distribution data, respectively. 
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log^oP [Tage] 

Fig. 15. The eccentricity-period diagram for the model of the pre-main sequence 
eigenevolution (A = 28, x ~ 0.75) at f = (upper panel) for stellar mass 0.1 < 
rrii/MQ < 1.1 and after cluster disintegration (bottom panel; note: Tage=days). 
Systems with semi- major axis after pre-main s equence eigenevolution a < 10 Rq 
have been merged. The upper envelope is from iDuguennov fc Mavorl l|l99ll ) - bi- 
naries are only observed to have e, IP below this envelope, the region above this 
envelope being forbidden, as pre-main sequence dissipational effects are understood 
to de-populate this region with 10^ yr. However, dynamical encounters in the clus- 
ter populate the eigenevolution region implying that short-lived forbidden orbital 
parameters should be observable in stellar clusters. Some of these are indicated as 
open circles. Eigenevolution (i.e. classical tidal circularisation) on the main sequence 
with Ams = 24.7, Xms ~ 8 applied to the data in the lower panel depopulates the 
eigenevolution region and circularises all orbits with period less than about 12 d. 
The dashed lines are constant peri-astron distances (eq. I168P for rperi = A Rq and 
rrisys = 2.2, 0.64, 0.2 Mq (in increasing thickness). Note that horizontal and vertical 
cuts through this diagramme produce eccentricity and period distribution functions 
(as well as mass-ratio distributions) for short and long-period systems in agreement 
with the observational constraints. Th e initial orbits with P > 10**'^ d come about 
due to crowding. From iKroupal l|l995dl ). 
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logjoP [days] 



Fig. 16. The period distribution functions (IPF: eq. [165] with ea. ll67l and for stellar 
masses 0.1 < mi/M© < 1.1). The dashed histogramme is derived from IPF using the 
eigenevolution and feeding algorithms, and represents the binary-population at an 
age of about 10"" yr, while the solid histogramme follows from the dashed one after 
passing through a cluster with initially A'^ — 200 binaries and ro.s — 0.8 pc. The 
agreement of the dashed histogramme with the observational pre-main sequence 
data (as in fig. I13p . and of the solid histogramme with the observational main- 
sequence (Galactic-field) data (also as in fig. I13|l is good, noting that the longest- 
period TTauri binary population is expected to show some disruption. A full model 
of the Galactic field late-type binary population has therewith been arrived at which 
unifies all available but apparently discordant observational data (see also figs 1141 
[T5]and[T7l). 
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Fig. 17. The mass-ratio distribution for stars witii 0.1 < m/M© < 1.1 is the 
solid histogramme, whereas the initial mass-ratio distribution (random pairing from 
the canonical IMF; after eigenevolution and feeding, at an age ~ 10^ yr) is shown 
as the dashed histogramme. The solid histogramme follows from the dashed one 
after passing through a cluster with initially TV — 200 binaries and ro.s = 0.8 pc. 
The observational data (solid dots, iReid fc GizisI 1 19971 ) have been obtained after 
removing WD companions and scaling to the model. This solar-neighb ourhood 8 pc 
samp le is not complete and may be biased towards g = 1 systems (jHenrv et al.l 
1 19971 ). Nevertheless, the agreement between model (solid histogramme) and the data 
is striking. A full model of the Galactic field binary population has therewith been 
arrived at which unifies all available but apparently discordant observational data 
(see also figs [Ml [T5l and [Te)) . 



Note that initial distributions are derived from birth distributions. This is 
to be understood in terms of these initial distributions being the initiahsa- 
tion of A^— body experiments, while the birth distributions are more related 
to the theoretical distribution of orbital parameters before dissipational and 
accretion processes have a major effect on them. The birth distributions are, 
however, mostly an algorithmic concept. Once the iV— body integration is fin- 
ished, e.g. when the cluster is dissolved, the remaining binaries can be evolved 
to the main-sequence distributions by applying the same eigenevolution algo- 
rithm above, but with parameters 



A„is = 24.7 , Xms = 8. 



(175) 
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The need for A„is < A and Xms > X to ensure for example the tid al 



circularisation period of 12 days for G dwarfs (jDuauennov fc Mavoij|l99l[ ). is 
nicely consistent qualitatively with the shrinking of pre-main-sequence stars 
and the emergence of radiative cores that essentially reduce the coupling be- 
tween the stellar surface, where the dissipational forces are most effective, and 
the center of the star. 



The reader is also directed to Mardling fc AarsethI ( 2001 ) where a model 



of tidal circularisation is introduced into the A''— body code. 

Finally, the above work and t he application to the ONC and Pleiades 



(jKroupa. Aarseth fc Hurlev II2001I ) suggests the following hypothesis: 



Initial Binary Universality Hypothesis: The initial period fea. ll67p . ec- 
centricity (eq. I155p , and mass-ratio (random pairing from canonical IMF) dis- 
tributions constitute the parent distribution of all late- type stellar populations. 



Can this hypothesis be disproven? 

4.3 The initial binary population massive stars 

The above semi-empirical distribution functions have been formulated for late- 
type stars (primary mass nii ^ 1 Mq) as it is for these that we have the best 
observational constraints. It is not clear yet if they are also applicable to 
massive binaries. 

An approach taken bv lClarke fc Prinde I (1199 2") is to constrain the binary 



properties of OB stars by assuming that runaway OB stars are ejected from 
star-forming regions. About 10-25 per cent of all O stars are runaway stars, 
while about 2 per cent of B stars are runaways. This approach leads to the 
result that massive stars must form in small— Af groups of binaries that are 
biased towards unit mass ratio. This is a potentially powerful approach, but 
it can only constrain the properties of OB binaries when they are ejected 
which occurs after substantial dynamical encounters in the cluster core which 
typically lead to the mass-ratio evolving towards unity as the involved binaries 
harden. The true birth properties of massive binaries therefore remain obscure, 
and one needs to resort to A^— body experiments to test various hypothesis 
given the observational constrains. One such hypothesis could be, for example, 
to assume massive stars form in binaries with birth pairing properties as for 
low- mass stars fsection l4.2l) . i.e. most massive primaries would have a low-mass 
companion, and to investigate if this hypothesis leads to the observed number 
of runaway massive stars through dynamical mass segregation to the cluster 
core and partner exchanges through dynamical encounters there between the 
massive stars. 

Apart from the fraction of runaway stars, direct surveys have lead to some 
insights as to the binary pro perties of the observed massive stars. Thus, for 
example. iBaines et al.l ( 2006 ) report a very high (/ « 0.7 ±0.1) binary fraction 



among Herbig Ae/Be stars with the binary fraction increasing with increasing 
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primary mass. Furthermore, they find that the circum-binary disks and the 
companions appear to be co-planar thereby supporting a fragmentation origin 
rather than collisions or capture as the origin of massive binaries. Most O stars 
are b elieved to exist as short-period binaries with g « 1 (jGarcfa fc Mermilliod 
20011 ). at least in rich clusters. On the other hand, small-g appear to be 
favoured in less substanti al clusters such as t he ONC, being consistent t here 



with random pairing (P reibisch et al. 19991 ). Kouwenhoven et al. ( 2005 ) re- 



port the A and late-type B binaries in the Scorpius 0B2 association to have 
a mass-ratio distribution not consis tent with random pa i ring. The lower limit 
on the binary fraction is 0.52, while [Kouwenhoven et aL (2007) update this to 
a binary fraction of 72 per cent. They also find that th e semi-major axis dis- 
tribut ion contains too many close pairs compared to a lDuquennov &: Mavoi 
( 199lh log- normal distribution. These are important constraints, but again, 
they are derived for binaries in an OB association, which is an expanded ver- 
sion of a dense star cluster (section II. 2[) and therefore hosts a dynamically 
evolved population. 

Given the above results, perhaps the massive binaries in the ONC represent 
the primordial population, whereas in rich clusters and in OB associations the 
population has already dynamically evolved through hardening and compan- 
ion exchanges to that observed there {fq rising towards q = 1). This possibility 
needs to be investigated using high-precision A'^-body computations of young 
star clusters. The first, simplest hypothesis to test would be to extend the 
pairing rules of section 14.21 to all stellar masses, perform many (because of 
the small number of massive stars) A'^— body renditions of the same basic pre- 
gas expulsion cluster, and to quantify the properties of the emerging stellar 



population at various dynamical times (jKroupa Il200lh . 



Another approach would be to constrain a and TO2 for a given mi ^ 5 Mq 
such that 

Sbin ~ (176) 

(eq. I14ip . Or one can test the initial massive-star population given by 



a < 



(177) 



which follows from stating that the density of a massive binary, 2 x 3/(a'^ 4 tt), 
be larger than the cluster-core density, A"ob 3/(r^ 4 tt) . 

So far, none of these possibilities have been tested, apart from extending 
the Initial Binary Universality Hypothesis (p. [TT]) to massive stars 



(|Kroum 1120011) . 
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5 Summary 

The above material gives an outline of how to set up an initial, birth or pri- 
mordial stellar population such that it resembles observed stellar populations. 
In section |4j2] a subtle differentiation was performed between initial and birth 
populations, in the sense that an initial population is derived from a birth 
population through initial processes that act too rapidly to be treated by an 
A'^— body integration. 

An A'^— body stellar system is generated for numerical experiments by spec- 
ifying its 3D structure and velocity field (section[21) , the mass distribution of its 
population (section [3]) and the properties of its binary population (section 2]). 

Given the distribution function discussed here, and the existing numerical 
results based on these, it is surprising how universal the stellar and binary 
population turns out to be at birth. A dependence of the IMF or the birth 
binary properties on the physical properties of star-forming clouds cannot be 
detected conclusively. In fact, the theoretical proposition that there should 
be a dependency is falsified, except possibly (i) in the extreme-tidal field 
environment at the Galactic centre, or (ii) in the extreme proto-stellar density 
environment of ultra-compact dwarf galaxies, or (iii) for extreme physical 
environments (pp. EIHMI)- 

The unified picture that has emerged concerning the origin of stellar pop- 
ulations is that stars form according to a universal IMF and mostly in binary 
systems and in very dense clusters that expel their residual gas and rapidly 
evolve to T- or OB-associations. If the latter are massive enough, the dense 
embedded clusters evolve to populous OB associations that may be expanding 
rapidly and that contain cluster remnants which may reach globular cluster 
masses and beyond in intense star-bursts. This unified picture naturally ex- 
plains the high infant weight loss and infant mortality of clusters, the binary 
properties of field stars, possibly thick disks of galaxies and the existence 
of population II stellar halos around galaxies that have old globular cluster 
systems. 

Some open questions that remain are quite obviously, why the star- 
formation product is so universal (within current constraints), and how mas- 
sive stars are distributed in binaries and if they form at the centres of their 
clusters, why the cluster mass of « 10^ Mq is special, and which star-cluster 
population is a full solution to the inverse dynamical population synthesis 
problem (p. [7T|) . Naturally, many more observations are required not only of 
topical high-redshift star-burst systems, but also of the more mundane low- 
redshift and preferably local star-forming objects and globular and open star 
clusters, to further refine the above broad picture. 
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